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ABSTRACT 

The  thesis  compares  the  analytical  solution,  tyvo  marine  classification  society  design 
rules,  and  two  design  guides  against  experimental  results  for  predicting  the  failure  modes 
(general  instability,  axisymmetric  buckling,  and  asymmetric  collapse  of  the  shell  [lobar 
buckling])  and  failure  pressures  of  ring-stiffened  cylinders 

The  analytical  solution  is  first  summarized  based  on  several  sources.  The  design  rules  for 
the  classification  societies  and  the  design  guidance  from  two  sources  are  then  presented  with 
brief  explanations  for  each  one.  The  design  rules  used  are:  American  Bureau  of  Shipping  {Rules 
for  Building  and  Classing  Underwater  Vehicles,  Systems,  and  Hyperbaric  Facilities,  1990)  and 
Germanischer  Lloyd  {Rules  for  Underwater  Technology,  1988).  The  design  guides  used  were 
Society  of  Naval  Architects  and  Marine  Engineers  {Submersible  Vehicle  Systems  Design,  1990) 
and  Massachusetts  Institute  of  Technology  Course  13A  Professional  Summer  Notes  (MIT  13A 
Submarine  Design  Trends,  2001). 

The  United  States  Navy  Naval  Sea  Systems  Command,  Submarine  Structural  Integrity 
Division  supplied  experimental  data  for  four  cylinders  that  covered  the  failure  modes  and 
allowed  comparison  between  experiment  and  design  rules  /  guidance. 

The  comparison  of  experimental  to  predicted  data  found  that  the  design  codes  and  design 
guides  performed  adequately  in  predicting  axisymmetric  yield  and  asymmetric  buckling.  The 
performance  of  the  design  codes  and  guides  in  predicting  failure  by  general  instability  was 
unsatisfactory.  For  the  experimental  failures  by  general  instability,  the  design  codes  and  guides 
predicted  significantly  higher  failure  pressures  than  those  experimentally  determined;  resulting  in 
the  design  codes  and  guides  actually  predicting  failure  by  axisymmetric  yield  in  stead  of  general 
instability.  These  inconsistencies  in  the  predictions  of  failure  mode  and  pressures  for  general 
instability  should  be  further  explored  to  determine  causes  and  corrections. 
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Chapter  1:  Introduction 


The  widespread  use  of  stiffened  cylinders  in  the  marine  industry  has  generated  many 
studies  into  the  stability  and  failure  of  these  cylinders  and  methodology  for  failure  prevention. 

Of  primary  concern  to  entities  involved  with  the  use  of  manned  submersible  vehicles  is  the 
design  of  ring-stiffened  cylinders;  this  type  of  stiffened  cylinder  is  used  for  significant  portions 
of  the  pressure  hull.  Over  a  hundred  years  of  theoretical  and  experimental  research  has  led  to  a 
general  understanding  of  the  mechanics  of  failure  for  these  cylinders.  Based  on  this  research, 
marine  Classification  Societies,  such  as  the  American  Bureau  of  Shipping  (ABS),  the  American 
Petroleum  Institute  (API),  NORSOK  and  Germanischer  Lloyd  (GL)  have  promulgated  design 
rules  to  provide  guidelines  on  the  design  and  building  of  stiffened  cylinders  for  safe  operation. 
Other  design  and  analysis  theories  and  guidance  are  available  in  texts  such  as  those  published  by 
the  Society  of  Naval  Architects  and  Naval  Engineers  (SNAME). 

1.1  Failure  Types 

There  are  three  primary  types  of  failure  of  ring-stiffened  cylinders.  They  are 
axisymmetric  yielding  (AY)  of  the  shell  between  stiffeners,  asymmetric  buckling  of  the  shell 
between  stiffeners  (Lobar),  and  general  instability  of  the  shell  and  stiffeners  (GI).  Axisymmetric 
yield  generally  occurs  when  the  shell  is  relatively  heavy  and  the  frames  are  closely  spaced. 

Lobar  buckling  can  occur  when  the  shell  is  relatively  light  and  the  frames  are  strong  and  widely 
spaced.  General  Instability  can  occur  when  the  cylinder  is  relatively  long,  the  shell  is  thin  and 
the  frames  are  small.  General  Instability  is  very  dependent  upon  eccentricities  in  the  shell,  which 
tend  to  lower  the  cylinder’s  resistance  to  the  General  Instability  mode.  [1] 
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As  analyzed  here  General  Instability  is  presumed  to  occur  in  the  elastic  region  of  the 
stress  -  strain  curve  of  the  material.  Cylinders  also  fail  by  inelastic  General  Instability,  which 
occurs  at  significantly  lower  pressures  than  that  of  elastic  General  Instability.  Failure  by  this 
mode  is  not  addressed  by  the  design  rules.  Other  modes  of  failure  also  exist  such  as  multi- wave 
instability,  which  is  a  sub-type  of  General  Instability.  It  can  occur  in  both  the  elastic  and 
inelastic  regions.  Again  the  design  rules  do  not  address  this  failure  mode.  Several  of  the 
classification  societies  address  stiffener  tripping,  which  is  the  rotation  of  a  stiffener  away  from 
perpendicular  with  the  shell,  however  stiffener  tripping  is  usually  a  precursor  to  general 
instability  and  is  not  a  separate  major  failure  mode. 

1.2  Concept  Exploration 

For  this  thesis,  an  emphasis  was  placed  on  exploring  how  the  various  design  rules 
predicted  failure  of  cylinders  that  replicate  modern  submarine  design  (i.e.  the  shell  was  relatively 
thick  compared  to  the  diameter  of  the  cylinder).  This  was  facilitated  by  experimental  failure  data 
from  the  U.S.  Navy’s  Naval  Sea  Systems  Command  (NAVSEA)  Submarine  Structural  Integrity 
Division.  The  analysis  was  limited  to  ring-stiffened  cylinders  to  remain  consistent  with  the 
primary  concern  of  the  thesis.  The  test  cylinders  that  were  chosen  had  failed  in  all  three  possible 
modes,  allowing  for  comparison  of  the  design  rules  in  all  modes  of  failure.  . 

1.3  Analysis  Techniques 

For  the  analyses,  the  design  rules  for  the  various  classification  societies  and  design 
guidance  were  programmed  into  MATHCADtm  for  consistency  of  approach,  ease  of  symbolic 
representation,  and  quickness  of  computations.  The  analytical  solution  was  also  programmed 
into  MATHCAD™  for  comparison.  The  experimental  data  from  the  four  test  cylinders  were 
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then  input  into  each  computer  code.  The  codes  gave  a  failure  pressure  for  each  type  of  failure. 
The  lowest  calculated  pressure  was  then  considered  the  failure  pressure  and  the  corresponding 
mode  was  designated  the  failure  mode.  The  failure  modes  and  pressures  were  then  compared 
against  the  experimental  results  with  emphasis  placed  on  agreement  between  codes  and 
experiment  on  failure  mode  (first  priority)  and  then  failure  pressure  (second  priority).  If  a 
different  failure  mode  was  predicted  than  that  experimentally  found  for  a  particular  cylinder,  an 
analysis  was  performed  for  agreement  between  predicted  failure  pressures  and  also  the  closeness 
of  the  failure  mode  pressures. 

1.4  Design  Rules  and  Guides  Examined 

There  were  two  classification  society  design  rules  examined;  The  American  Bureau  of 
Shipping  (rules  from  Rules  for  Building  and  Classing  Underwater  Vehicles,  Systems  and 
Hyperbaric  Facilities,  1 990  Edition)  (reference  2)  and  Germanischer  Lloyd  {Rules  for 
Underwater  Technology,  1988  Edition)  (reference  3).  Design  guides  included:  The  Society  of 
Naval  Architects  and  Marine  Engineers  (SNAME)  {Submersible  Vehicle  Systems  Design,  1990) 
(reference  4)  and  the  Massachusetts  Institute  of  Technology  (MIT)  Course  13A  Professional 
Summer  (13 A  PS)  Submarine  Design  Trends  Course  notes  (reference  5)  was  also  used.  The  13A 
PS  is  based  upon  the  1967  version  of  the  SNAME  publication  Principles  of  Naval  Architecture 
(PNA)  with  some  modifications  determined  by  Harry  Jackson  (CAPT,  USN  Ret.).  A  third 
classification  society,  the  American  Petroleum  Institute  (API)  was  planned  to  be  used,  however 
the  immediately  available  rules  from  API  were  not  valid  for  the  experimental  data  used. 
However,  API  does  have  rules  to  cover  the  types  of  cylinders  examined  here. 

Analytic  solutions  for  the  three  failure  modes  were  gathered  from  several  sources.  These 
sources  included  Hydrostatically  Loaded  Structures:  The  Structural  Mechanics,  Analysis  and 
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Design  of  Powered  Suhmersihles  by  William  Nash  (reference  6),  Principles  of  Naval 
Architecture,  1967  edition  (reference  1),  a  David  Taylor  Model  Basin  technical  paper  by  J.G. 
Pulos  and  V.L.  Salerno  (reference  7),  and  several  journal  articles  from  the  Transactions  of  The 
Royal  Institution  of  Naval  Architects  by  S.  Kendrick  (reference  8)  and  C.  T.  F.  Ross  (reference 

9). 

In  Chapter  2  a  brief  discussion  on  the  terminology  of  ring-stiffened  cylindrical  shells  is 
produced  along  with  a  short  derivation  of  the  basic  stresses  in  these  cylinders.  Chapter  3 
contains  the  summary  of  analytic  solutions  used  to  predict  failure  (the  solutions  are  not  re¬ 
derived),  Chapter  4  has  the  summary  of  the  design  rules  and  guides  from  the  various  sources. 
Chapter  5  describes  the  test  cylinders  and  then  compares  the  results  from  experiment  to  the 
results  predicted  by  the  design  rules  and  guides.  Chapter  6  summarizes  the  results  with 
recommendations  for  further  study. 


Chapter  2:  Overview  of  Ring-Stiffened  Cylindrical  Shells 


Ring-stiffened  cylinders  are  the  prevalent  construction  type  in  the  mid-bodies  of  modem 
submersibles.  The  stiffeners  provide  additional  strength  to  the  shell  that  is  required  for  the 
pressure  differential  between  the  external  hydrostatic  pressure  and  the  internal,  approximate 
atmospheric  pressure.  [6] 

2.1  Terminology 

The  various  classification  societies  use  slightly  different  terminology  for  the  cylinder 
geometries  and  properties.  In  the  computerized  design  rules,  the  symbols  used  by  each 
classification  society  are  generally  used  to  avoid  confusion  between  the  published  code  and  the 
programs.  All  of  the  codes  require  computation  of  the  moment  of  inertia  (I)  of  the  ring  stiffener 
and  the  moment  of  inertia  (le)  of  an  effective  ring  stiffener  (the  frame  and  some  length  of 
attached  shell).  The  formulas  for  I  and  came  from  [10]. 

The  terminology  can  be  divided  into  two  categories:  material  properties  and  geometry. 
All  stresses  and  pressures  are  in  pounds  per  square  inch  (psi),  lengths  are  in  inches  (in),  areas  are 
in  square  inches  (in^)  and  moments  of  inertia  are  in  inches  to  the  fourth  power  (in'*). 

2.1.1  Material  Properties 

1)  Modulus  of  Elasticity  (E):  The  slope  of  the  linear  region  of  the  stress-strain  curve  of  a  given 
material.  For  all  of  the  calculations  a  value  of  3  X  10^  psi  (a  common  value  for  steel)  was  used. 

2)  Poisson’s  Ratio  (v  or  p):  The  ratio  of  lateral  strain  to  axial  strain  in  a  material.  For  all 
calculations,  a  value  of  0.3  (a  general  figure  for  steel)  was  used. 
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3)  Yield  Strength  (Oy  or  fy  or  k):  An  arbitrary  value  for  materials  marking  the  onset  of  plastic 
deformation  of  a  material.  This  is  usually  considered  the  point  of  0.2%  permanent  strain.  This 
value  is  of  primary  concern  in  the  axisymmetric  yield  calculations.  This  parameter  ranged  in 
value  from  65,500  psi  to  1 57,000  psi. 

2.1.2  Geometry 

2.1.2.1  Cylinder  Geometry  (see  Figure  1) 

1)  Cylinder  Length  (L  or  Lb  or  Lc):  Overall  length  of  the  cylinder  between  supports.  This  varies 
by  cylinder  and  is  of  primary  concern  for  General  Instability  calculations 

2)  Shell  Mid-plane  Radius  (R);  Radius  from  centerline  of  cylinder  to  the  shell  mid-plane. 

3)  Shell  plate  thickness  (t  or  tp  or  s  or  h):  Thickness  of  the  cylinder  shell  plating. 

2.1.2.2  Stiffener  Geometry  (see  Figure  1) 

1)  Length  between  Stiffeners  (L  or  Lf  or  Ls):  Distance  between  centerlines  of  adjacent  stiffeners. 
This  distance  was  assumed  to  be  constant  for  each  cylinder.  This  dimension  is  important  for 
both  Axisymmetric  Yield  and  Lobar  Buckling. 

2)  Web  Height  (hw):  Length  of  the  web  for  ‘T’  stiffeners  or  the  height  of  the  stiffener  for 
rectangular  stiffeners. 

3)  Web  Thickness  (U):  Thickness  of  the  web  for ‘T’  stiffeners  or  the  thickness  of  the  stiffener  for 
rectangular  stiffeners. 

4)  Flange  Breadth  (bf):  Width  of  the  flange  for  ‘T’  stiffeners.  This  value  is  set  to  zero  for 
rectangular  stiffeners. 

5)  Flange  Depth  (d):  Thickness  of  the  flange  for ‘T’  stiffeners.  This  value  is  set  to  zero  for 
rectangular  stiffeners. 
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6)  Stiffener  Depth  (H);  The  distance  from  the  shell  to  the  end  of  the  stiffener.  This  is  the 
stiffener  height  for  reetangular  stiffeners  and  the  web  height  plus  the  flange  depth  for  ‘T’ 
stiffeners. 

7)  Faying  Length  (b):  The  distance  of  contact  of  the  stiffener  to  the  shell.  It  is  equal  to  the  flange 
width  for ‘T’  and  rectangular  stiffeners;  if  T’-beam  or  wide-beam  (WF)  stiffeners  are  used,  then 
equal  to  the  flange  breadth  in  contact  with  the  shell. 

8)  Effective  Length  of  the  Shell  (Le):  A  length  of  shell  to  be  considered  as  part  of  the  combined 
stiffener  and  shell  used  in  the  General  Instability  analyses.  This  length  is  usually  a  fraction  of 
the  stiffener  spacing,  but  many  times  will  be  equal  to  the  stiffener  spacing. 

9)  Area  of  the  stiffener  (A):  Cross-sectional  area  of  the  stiffener. 

10)  Effective  Area  of  the  stiffener  (Aeff):  Cross-sectional  area  of  the  combined  stiffener  and 
effective  length  of  shell  or  a  modified  area  of  the  stiffener  (A)  based  on  the  location  of  the 
stiffener. 

1 1)  Moment  of  Inertia  (I):  Area  moment  of  inertia  of  the  dedicated  stiffener,  used  by  some  codes 
as  part  of  the  calculation  for  the  effective  moment  of  inertia. 

12)  Effective  Moment  of  Inertia  (Lff):  Area  moment  of  inertia  for  the  combined  stiffener  and 
shell,  used  for  General  Instability  calculations. 
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Figure  1:  Cylinder  and  Stiffener  Geometry 


2.2  Classification  of  Stiffeners 

There  are  several  methods  to  classify  stiffeners.  Of  importance  to  the  analysis  of 
cylinders  is  the  location  of  the  stiffener:  either  internal  or  external  to  the  cylinder  shell.  Both 
types  are  used  for  cylinder  construction;  slight  modifications  to  the  design  rules  (concerning 
effective  areas  and  moments  of  inertia)  are  required  based  on  the  location.  External  stiffeners  of 
equal  size  to  internal  stiffeners  require  more  material  because  of  the  greater  circumferential 
length  but  offer  the  advantage  of  freeing  spacing  within  the  cylinder  for  equipment  /  living  space. 

[5] 
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Another  classification  of  stiffeners  are  field  (light  /  non-heavy)  and  heavy  (king  frame) 
stiffeners.  Heavy  stiffeners  are  substantially  larger  than  the  field  stiffeners  and  are  used  to 
reduce  the  effective  length  of  the  cylinder  for  general  instability  concerns.  Field  stiffeners  are  of 
uniform  size,  shape,  and  spacing.  For  the  current  analyses  no  heavy  stiffeners  were  used. 

2.3  Stresses  in  Cylinders 

A  brief  discussion  of  stresses  in  cylindrical  shells  is  required  to  set  up  the  derivation  of 
the  structural  mechanics  in  support  of  the  analytical  solution.  As  a  starting  point,  a  cylinder  can 
be  considered  a  thin-walled  structure  (shell)  if  the  ratio  of  the  mid-plane  radius  to  the  wall 
thickness  is  greater  than  ten.  This  assumption  allows  the  determination  of  the  stresses  by  statics 
alone.[10]  All  of  the  cylinders  considered  are  treated  as  shells.  A  second  assumption  in  the 
analysis  is  to  consider  the  hydrostatic  pressure  as  constant  across  the  cylinder. 

Cylindrical  shells,  exposed  to  hydrostatic  pressure,  have  two  basic  stresses  imposed  by 
the  pressure:  hoop  stress  and  axial  stress.  The  equations  for  the  stress  are: 


1)  Hoop  (circumferential)  Stress: 

(1) 

2)  Axial  Stress: 

rr  -P^ 

CT  ' — 

"  It 

(2) 

Where  p  is  the  external  pressure,  R  is  the  shell  mean  radius  and  t  is  the  shell  thickness.  Figure  2 
shows  the  derivation  of  these  two  equations. 

The  addition  of  ring  stiffeners  to  the  base  shell  complicates  the  hoop  stress  analysis  by 
introducing  non-uniform  deformation  of  the  shell  in  the  radial  direction.  There  is  also  a  beam- 
column  effect  due  to  the  pressure  acting  in  the  axial  direction.  These  effects  will  be  addressed  in 
Chapter  3. 
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Figure  2:  Basic  Cylindrical  Shell  Stresses 


1)  Axial  Stress 

Axial  force  on  a  vertical  section  of  the  cylinder  must  equal  the  axial  stress  times  the 
circumferential  area. 

Force  balance:  p-ti  R^  =  2  n-R-t 
Which  results  in: 

PR 

^axial  “  2.x 

2)  Circumferential  (Hoop)  Stress 

The  force  on  a  transverse  section  of  shell,  with  width  b  must  equal  the  hoop  stress  times  the  shell 
area. 

Force  balance:  p-2  R  b 
Which  results  in: 

PR 

®^hoop  •“  j 
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Chapter  3:  Analytic  Solutions 


Analytic  solutions  have  been  proposed  and  proven  for  the  three  major  failure  modes  of 
stiffened  cylinders.  This  chapter  briefly  describes  the  theories  and  equations  used  in  the 
experimental  analyses. 

3.1  History  of  Analyses 

The  failure  of  cylinders  exposed  to  external  pressure  has  been  studied  for  almost  1 50 
years.  The  first  attempts  at  understanding  cylinder  behavior  was  done  by  experiment  and 
empirical  relationships  in  the  1850’s.[6]  The  first  analytic  solution  for  a  non-reinforced  cylinder 
was  presented  by  G.  H.  Bryan  in  1 888.  [9]  The  first  analysis  of  a  reinforced  cylinder  appeared  in 
1913  by  R.  V.  Southwell,  followed  a  year  later  by  a  solution  to  the  elastic  buckling  of  a  thin  shell 
proposed  by  von  Mises.[6]  In  1934  Widenburg  proposed  a  solution  for  asymmetric  buckling  that 
is  independent  of  the  number  of  lobes  of  failure.  [6]  This  equation  is  the  one  used  in  the  current 
analysis.  Solutions  for  axisymmetric  yield  were  first  put  forward  by  von  Sanden  and  Gunther  in 
1920. [9]  Viterbo  presented  a  modified  version  of  Sanden  and  Gunther’s  solution  in  1930. [9] 
Pulos  and  Salerno  presented  a  solution  that  included  the  Sanden  and  Gunther  solution,  the 
Viterbo  modification  and  a  term  to  account  for  the  bending  stress  in  the  cylinder  caused  by  the 
axial  pressure.  [9]  The  Pulos  and  Salerno  solution  is  used  in  this  thesis.  For  elastic  general 
instability,  the  first  reported  analysis  was  presented  by  Tokugawa  in  1929.  In  1954  A.  R.  Bryant 
developed  a  similar  equation  using  a  different  methodology.  [6] 

Analytical  work  from  the  1950’s  onward  has  focused  on  obtaining  solutions  for  different 
boundary  conditions  and  more  fully  reconciling  the  analytic  predictions  with  experimental 
results  and  to  more  fully  understand  the  effects  of  initial  imperfections  in  the  cylinder’s  material 
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and  geometric  properties.  With  the  advent  of  power  digital  computers  and  the  use  of  finite 
element  analysis  there  has  been  great  strides  made  in  understanding  the  failure  of  cylinders. 

3.2  Current  Theory 

For  this  thesis,  analytic  solutions  were  collected  from  several  sources  for  the  modes  of 
failure.  A  comprehensive  theoretical  solution  that  addresses  all  modes  of  failure  is  not  presented. 
Reference  6  provides  a  good  summation  of  the  currently  used  analytic  solutions. 

A  first  indicator  of  the  failure  mode  of  a  cylinder  is  found  by  plotting  the  cylinder’s 
slenderness  ratio  (X)  against  the  ratio  of  the  shell  buckling  pressure(pc)  to  the  hoop  pressure  at 
yield(py)  (\i/).[l]  The  equations  for  these  factors  are  given  below. 


1 


Figure  3  shows  the  plot  of  vj/  versus  X.  If  the  slenderness  ratio  is  less  than  approximately  1.14  the 
cylinder  should  fail  by  axisymmetric  yield;  if  is  greater  than  1.14  then  the  cylinder  should  fail 
by  lobar  buckling.  If  the  shell  and  stiffeners  are  not  of  sufficient  size,  the  cylinder  may  fail  by 
general  instability  at  a  pressure  less  than  that  found  in  Figure  3.[1]  By  using  two  assumptions 
(the  material  is  steel  with  v  =  0.3  and  Lj/D  »  tp/D)  it  can  be  shown  that:  [1] 

1.30 

'"“-F  <5) 

This  is  the  buckling  part  of  the  curve  in  Figure  3. 
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Figure  3:  Failure  Pressure  Ratio  versus  Slenderness  Ratio 


A 


The  treatment  of  boundary  /  end  conditions  of  the  cylinder  is  a  vitally  important  factor  in 
the  analytic  solutions.  The  literature  is  full  of  discussion  on  what  types  of  end  condition  to  use, 
with  the  choices  ranging  from  full  clamped  ends  to  simply  supported  ends.  The  extreme  cases 
are  hard  to  create  in  reality  and  therefore  the  experimental  results  tend  to  fall  between  the  ranges 
of  predictions.  Experiments  have  shown  that  partially  clamped  cylinders  provide  significantly 
higher  failure  pressures  than  that  predicted  by  mathematical  models  utilizing  simply  supported 
ends.[8]  For  this  thesis,  no  discrete  boundary  conditions  were  required  to  be  stated  for  input  into 
the  equations. 

3.2.1  Axisymmetric  Yield 

As  mentioned  in  section  3.1  axisymmetric  yield  has  been  studied  since  the  1920’s.  The 
solution  summarized  here  was  put  forward  by  Pulos  and  Salerno  in  1961.  It  is  based  on  the 
previous  works  of  van  Sunden  and  Gunther  and  Viterbo  and  includes  a  previously  not  included 
“beam  -  column”  effect  due  to  the  hydrostatic  pressure  acting  in  the  axial  direction  of  the 
cylinder. [7]  The  governing  differential  equation  is: 
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(6) 


+ 


P-R 

2  dx^ 


Et 

w  +  — 


R 


2 


w 


Where: 


E  t 

0  — -  Y)  jg  flexural  rigidity  of  the  shell 

12(1  -v^J 

The  term  represents  the  beam  -  column  effect.  It  makes  the  solution  to  equation  (6)  a  non¬ 
linear  function  of  pressure  and  was  the  term  neglected  in  the  earlier  analyses  of  axisymmetric 
yield.  Figure  4  shows  the  coordinate  system  used  in  reference  7  to  derive  the  governing 
equation,  x,  (p,  and  r  are  the  axial,  circumferential  and  axial  coordinates  respectfully  with  u,  v, 
and  w  being  the  corresponding  displacements.  [7] 


Figure  4:  Coordinate  System  for  a  Cylindrical  Shell 
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Following  typical  practice  in  the  solution  of  non-homogeneous  differential  equations  the 
general  solution  of  the  governing  equation  can  be  written  as  the  sum  of  the  solution  of  the 
homogeneous  equation  and  a  particular  solution.  [6]  The  solution  to  the  homogeneous  equation 
produces  four  roots  (X,i-X4).  By  analysis,  placement  of  the  origin  of  the  coordinate  system  to  take 
advantage  of  symmetry,  and  trigonometric  identities,  the  general  solution  can  be  given  as: 

w=  5cosh/l,x  +  Fcosh/l,x-^^^(l-— )  (7) 

Et  2 

where  B  and  F  are  arbitrary  constants  which  can  be  found  by  applying  boundary  conditions  to 
the  equation.  [7].  After  further  mathematical  substitutions  several  dimensionless  parameters  were 
introduced  into  the  solution  to  allow  easier  solving  of  the  problem.  Four  of  these  parameters  (Fi- 
F4)  were  transcendental  functions  based  on  the  geometry  of  the  cylinder.  These  functions  were 
originally  graphed  to  allow  for  a  relatively  quick  solution  to  be  found  for  a  given  cylinder. 

Finally  an  equation  for  the  failure  pressure  of  a  given  cylinder  was  determined.  This  equation 
along  with  the  dimensionless  parameters  is  given  below. 
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A  ~ 


a  +  p  +  (l-p).F, 


An  iterative  process  is  required  for  the  general  case  where  the  parameter  y  is  not  zero.  The 
process  begins  with  assuming  that  y  is  zero,  finding  the  failure  pressure  then  recalculating  y  and 
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solving  the  equations  for  the  failure  pressure  again.  Usually  only  two  to  three  iterations  are 
needed  for  satisfactory  convergence  of  the  failure  pressure.  [7] 

3.2.2  Asymmetric  Buckling  (Lobar  Buckling) 

Asymmetric  buckling  is  characterized  by  circumferential  lobes  between  ring  stiffeners. 
As  noted  above  this  mode  of  failure  will  occur  when  the  cylinder’s  slenderness  ratio  is  relatively 
high.  This  can  be  further  characterized  by  a  relatively  thin  shell  thickness  and  widely  spaced 
stiffeners.[l]  R.  von  Mises  first  proposed  a  solution  to  the  buckling  of  un-stiffened  cylinders 
under  hydrostatic  pressure  in  1929.  He  assumed  sinusoidal  displacements  in  the  axial  and 
circumferential  directions  to  enable  solving  of  a  set  of  linearized  partial  differential  equations 
that  represented  the  elastic  action  of  the  shell.  [6]  He  eventually  obtained  the  following  equation 
for  the  buckling  pressure: 


The  buckling  pressure  is  dependent  upon  the  number  of  circumferential  lobes  (n),  which  must  be 
an  integer  value.  This  fact  requires  an  iterative  process  of  varying  n  until  the  lowest  pressure  is 
determined. 

In  1933  Widenburg  solved  the  above  equation  in  way  that  was  independent  of  n.  From 
test  data  the  buckling  pressures  from  the  von  Mises  equation  and  the  Widenburg  approximation 
differ  by  no  more  than  3.5%.  [6].  Further  investigation  into  the  bucking  of  stiffened  cylinders 
determined  that  the  Widenburg  equation  worked  very  well  by  replacing  the  length  of  the  cylinder 
by  the  length  between  stiffeners.  Therefore  the  Widenburg  equation  (equation  10),  shown 
below,  is  the  equation  used  in  the  current  analyses. 


3.23  General  Instability 

General  instability  consists  of  the  yielding  of  both  the  cylindrical  shell  and  ring  stiffeners. 
A  cylinder  may  be  susceptible  when  the  stiffeners  are  undersized  when  compared  to  the  shell 
thickness  and  the  cylinder  is  relatively  long. [6]  General  instability  may  occur  in  either  the  elastic 
or  inelastic  stress  region  of  a  material.  Elastic  general  instability  is  the  mode  covered  by  the 
available  literature  and  is  addressed  here.  Inelastic  general  instability  has  been  studied  mainly  by 
government  laboratories  and  most  of  the  knowledge  of  this  failure  mode  is  classified  material 
and  therefore  unavailable  to  the  present  author.  [11] 

The  first  analysis  of  general  instability  was  presented  by  T.  Tokugawa  in  1929.  [6]  His 
methodology  considered  the  failure  as  a  combination  of  the  failure  of  the  ring  stiffeners  and  shell 
buckling  with  each  taking  place  separately.  [9]  In  the  1940’s  S.  Kendrick  used  a  strain  energy 
methodology  to  determine  the  failure  pressure.  In  1 953  A.  R.  Bryant  used  a  simpler  strain 
energy  method  and  determined  an  equation  similar  to  both  Kendrick’s  and  Tokugawa’ s.  [1,9] 
The  Bryant  equation  is  therefore  used  for  this  analysis  and  is  shown  below  as  equation  10. 
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The  first  term  can  be  considered  the  failure  of  the  shell  and  the  second  term  can  be 
considered  the  failure  of  a  combined  stiffener  and  an  effective  length  of  shell.  [6]  This  effective 
length  has  had  much  discussion  in  the  literature  over  the  decades.  Bryant  assumed  the  length  to 
be  the  spacing  between  the  stiffeners,  but  others  have  proposed  various  corrections  based  on  the 
cylinder’s  geometric  and  material  properties. [5]  For  this  thesis,  the  effective  length  used  was  the 
stiffener  spacing. 

Similar  to  the  von  Mises  buckling  pressure  determination,  the  number  of  circumferential 
lobes  (n)  must  be  found  that  minimizes  the  failure  pressure.  The  number  of  circumferential  lobes 
is  usually  between  2  and  4.  [5]  The  factor  A,  is  the  number  of  longitudinal  lobes  in  the  cylinder. 
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Chapter  4:  The  Design  Rules  and  Analysis  Tools 


The  two  classification  societies’  design  rules  that  were  utilized  were  the  American 
Bureau  of  Shipping  and  the  Germanischer  Lloyd  rules.  These  were  chosen  for  their  availability, 
their  different  levels  of  simplification  of  equations,  their  coverage  of  the  specific  geometries  of 
the  experimental  cylinders,  and  their  inclusion  of  all  three  failure  modes.  The  Society  of  Naval 
Architects  and  Marine  Engineers  Submersible  Vehicle  Systems  Design  and  Principles  of  Naval 
Architecture  (PNA)  were  chosen  for  analysis  as  SNAME  is  the  primary  design  society  in  the 
United  States  and  has  comprehensive  guidelines  for  cylinders.  The  MIT  13A  Professional 
Summer  Submarine  Design  Trends  notes  were  used  because  of  the  author’s  familiarity  and  the 
complete  analysis  of  the  failure  modes. 

4.1  American  Bureau  of  Shipping  (ABS) 

The  ABS  design  rules,  as  delineated  in  the  Rules  for  Building  and  Classing  Underwater 
Vehicles,  Systems  and  Hyperbaric  Facilities  (1990),  give  a  brief  and  conservative  approach  for 
determining  the  critical  /  collapse  pressures  for  each  failure  mode.  The  ABS  design  rules  do  not 
explicitly  name  the  failure  modes,  but  there  are  distinct  equations  for  the  three  modes. 

4.1.1  Axisymmetric  Yield 

This  mode  is  designated  the  yield  pressure  at  midbay  and  midplane  of  a  cylinder.  [2]  The 
formula  accounts  for  the  major  parts  of  the  analytic  solution,  but  uses  single  value  functions  for 
the  shell  parameters  and  does  not  explain  the  functions  of  each  part  of  the  equation.  The  rules  do 
account  for  the  difference  between  internal  and  external  frames  by  squaring  the  mid-plane  radius 
to  stiffener  radius  ratio  for  external  stiffeners.  The  base  equation  follows: 


27 


(12) 


t 


Where 


F:= 


A  +  t^t  + 


2NtL 

e 


Internal  Stiffeners  or 


2 

External  Stiffeners 


N  —  ~  cos(2  Q) 

sinh(2Q)  +  sin(2Q) 

^  (sinh(Q)  cos(Q)  +  cosh(Q)  sin(Q)) 
sinh(2Q)  +  sin(2Q) 

sinh(2Q)  -  sin(2  Q) 
sinh(2Q)  +  sin(2Q) 


4.1.2  Asymmetric  Buckling 

This  is  called  the  Von  Mises  buckling  pressure  for  a  cylinder.  It  is  the  Widenburg 
approximation  that  is  used  by  most  of  the  classification  societies. 
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(l-v^) 


2.42E- 


4  ,  /  \2 

• - 0.45f— ] 

2-R  I  2-R  J 


The  ABS  code  then  has  a  range  of  allowable  pressures  depending  upon  the  ratio  of  Pm  to 
Py.  Below  is  the  logic  for  the  maximum  allowable  working  pressure  for  inter-stiffener  strength 
(Pc).  Further  safety  factors  are  then  applied  to  lower  the  allowable  pressure  in  practice. 


p  P 

m 

—  if  —  <  I 
2  P 

y 

p  ^  p 

y  m 

P^-  1 - —  if  1  <  —  <3 

y  2p  p 


5  „  ..  m  . 
-P  if >3 

6  y  Py 


4.1.3  General  Instability 

For  the  elastic  general  instability,  the  ABS  code  uses  the  Bryant  equation  with  an 
effective  length  not  equal  to  one  frame  space.  The  effective  length  is  given  below. 


L  :=  mil/  *  ^  The  equation  for  the  failure  pressure  is  broken  into  three  parts,  but  the 

'  I 


total  is  equivalent  to  the  Bryant  formula. 


Et  EvIgA2(n) 

Pni(")  :=  ^  — ; — 

^  R^L 


Where; 
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4.2  Germanischer  Lloyd 

The  design  rules  for  Germanischer  Lloyd  are  from  Rules  for  Underwater  Technology 
(1988  edition).  The  rules  address  all  three  failure  modes  very  thoroughly  and  flow  charts  are 
provided  to  aide  in  programming  the  code  for  computer  use.  These  rules  also  address  out-of- 
roundness  up  to  a  nominal  value  of  0,5%  of  the  mid-plane  shell  radius  by  determining  a  reduced 
allowable  pressure. [3] 


4.2.1  Axisymmetric  Yield 

The  code  for  axisymmetric  yield  resembles  the  Pulos  and  Salerno  methodology  with 
some  additions  to  account  for  the  transition  of  the  material  into  the  plastic  range  during  yield. 
The  methodology  consists  of  first  guessing  a  pressure  lower  than  the  failure  pressure  (the  test 
cases  used  1  psi  as  the  starting  point)  then  iterating  through  a  series  of  equations  to  determine  the 
shell  stresses.  If  the  determined  shell  stress  becomes  greater  than  0.8  of  the  yield  strength  of  the 
material,  then  a  Secant  Modulus  of  Elasticity  (Es),  a  Tangent  Modulus  of  Elasticity  (E.)  and  a 
plastic  range  Poisson  s  Ratio  (Up )  is  calculated  and  used  to  determine  the  other  calculation 
factors.  The  determination  of  the  integer  m  in  equation  (14)  is  described  in  the  code  below. 


Where: 


Must  iterate  on  the  integer  m  until  trial  is  <  trial2.  By  meeting  this  condition  the 
minimum  failure  pressure  is  found  (with  the  minimum  m). 

ajLi 

trial  := - 

n 

triaI2  := 

H,:=[i  +  H4{h2^-3(i-Vp^)]] 

H2  [(2 ^-pK] 

H3:=[(l-2.Vp)-(2-Vp).Kj 

,1-  ^ 


2 

H2H4 

»1  j 
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H3H4 


C3  *—  I  + 


H2H3H4 


C4  is  not  used  for  cylinders  without  heavy  stiffeners. 


C5 a-Lj 


Past  • 


2-s^E 

-  v^) 


critical  pressure,  elastic,  calculation  factor 


p  =  1  psi  (arbitrary  low  pressure) 

This  is  the  actual  starting  point  of  the  iteration.  A  pressure  (p)  is  guessed;  the  code  then 
calculates  the  failure  pressure.  If  the  failure  pressure  is  greater  than  the  guessed  pressure,  the 
applied  pressure  is  increased  by  a  set  amount.  This  continues  until  the  two  pressures,  p  and  pm 
are  within  a  chosen  delta  of  each  other.  This  is  essentially  the  method  of  Pulos  and  Salerno  with 
a  simpler  method  for  iterating  on  the  presumed  pressure. 


C.  :=  - 


+  G 


Cg Cj  Cg 


2  /  s-L, 


A  b  f  b  ) 

H  I-  1 - Fi 

s-L,  Lj  1  L,  I 
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Cji:- 


0.91 


1  -  V 


2  2 
cosh  -  cos 


^  cosh^Cg)*sinh(Cg)  cos^C^^-sin^C^^ 


cosh  ^Cg)'Sin(Cp)  sinh^Cg)-cos 


cosh  (Cg^'Sinh^Cg)  cos  ♦sin 


cos  sin  cosh  ^Cg^  sinh(Cg) 


^  2  cosh(Cg)-sinh(Cg)  cos  (^)  •sin^C^) 


cosh  ^Cg^  *  sin^C^^  sinh  ^Cgj*  cos 


F4:= 


2  cosh ^Cg)- sinh ^Cg^  cos(c^'S\n(c^ 

1  ”  ^  .  . . . . - _  .1.  ,  , 


The  functions  FI  through  F4  are  equivalent  to  those  presented  in  the  analytic  solution. 


a  := 


shape  factor 


The  following  set  of  equations  determines  the  stresses  at  the  mid-bay  of  the  shell. 


•= 


-pR 


2  ^10' 


'^xT-V  2"^10'^lf*^4 
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If  the  calculated  stress  is  >  0.8k  then  the  following  equations  are  used  to  determine  the 
material  properties  in  the  elastic  -  plastic  region.  The  code  is  somewhat  circular,  the  assumption 
was  made  to  determine  the  strain  (ei)  from  the  previously  calculated  stress  then  to  find  a  new 
stress  level  by  the  equation  given  in  the  code. 
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4.2.2  Asymmetric  Buckling 


For  asymmetric  bucking  the  Germanischer  Lloyd  code  uses  a  modified  version  of 
equation  9.  Therefore  the  dependence  upon  the  number  of  lobes  is  not  removed  and  should  give 
a  better  estimate  of  the  buckling  pressure.  The  equation  is  given  below. 


Pnl(n) : 


EsPjjj(n) 


Where: 


Pnl(n):=- 


s  -  1  +  A,j 

-  +  1  +  - 7 - ^ 


2  ,  ..2 
n  —  1  +  .5'A.| 


4.2.3  General  Instability 

For  elastic  general  instability,  the  code  uses  the  Bryant  equation  with  the  effective 
moment  of  inertia  based  upon  the  effective  length  of  shell  as  defined: 

^test  ■=  b  +  yj2  R-s 

Lg  :=  Lgjgjj  if  Lgjgjj  <  Lj  where  Li  is  the  frame  spacing 
Lj  otherwise 

The  base  equation  for  general  instability  also  uses  a  modified  radius  (Ro)  which  is  the  radius  to 
the  centroid  of  the  combined  stiffener  and  effective  length  of  shell.  For  this  thesis  no  “heavy”  or 
“king  frame”  stiffeners  were  used,  however  the  code  allows  for  these  stiffeners  with  several 
more  equations. 


Pnla(")  ■=Po(")  +  Pl(") 
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Where: 


4.3  Society  of  Naval  Architects  and  Marine  Engineers  (SNAME) 

The  SNAME  code  was  taken  from  Submersible  Vehicle  Systems  Design  (1990).  SNAME 
is  not  a  classification  society  as  ABS  and  Germanischer  Lloyd,  but  it  does  provide  many  of  the 
professional  and  technical  resource  material  for  naval  architects  in  the  United  States.  The 
current  codes  are  revised  versions  of  those  found  in  the  SNAME  publication  Principles  of  Naval 
Architecture  (1967  edition),  reference  1 .  The  code  addresses  all  three  modes  of  failure  with  the 
least  complex  set  of  equations  of  those  under  current  study. 

4.3.1  Axisymmetric  Yield 

For  axisymmetric  yield,  the  methodology  simplifies  the  analytic  solution  by  using  the 
older  theory  of  von  Sanden  and  Gunther  instead  of  the  solution  put  forward  by  Pulos  and  Salerno 
in  reference  7.  The  methodology  thus  neglects  the  beam-column  effect  of  the  axial  pressure. 

This  is  justified  by  designing  the  shell  to  yield  vice  buckle.  [1] 
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4.3.2  Asymmetric  Buckling 

For  Lobar  buckling,  the  code  uses  the  Widenburg  approximation  equation  that  is 
independent  of  the  number  of  failure  lobes.  This  is  the  same  as  used  by  ABS. 


37 


4.3.3  General  Instability 


For  general  instability,  the  code  uses  the  unmodified  Bryant  equation;  therefore  the 
length  of  shell  used  in  the  moment  of  inertia  calculation  is  one  stiffener  spacing.  The  code 
suggests  that  the  critical  value  of  n  is  between  2  and  4. 


Pcri(") 


R 


m 


(n^  -  i)  e  L 


2  ,  m 

n  -  1  +  — 
2 


(  2  l] 
\n  +  m  / 


R  Lr 


(20) 


Where: 


m:= 


7cR 

h 


4.4  MIT  13A  Professional  Summer  Submarine  Design  Trends 

The  Professional  Summer  notes  are  a  compilation  of  design  theories  from  Harry  Jackson 
(CAPT,  USN,  Ret.)  used  to  instruct  the  MIT  course  13A  students  in  the  basics  of  submarine 
design.  The  actual  purpose  of  the  code  is  to  evaluate  a  proposed  submarine  design  given  depth 
and  material  criteria;  some  modifications  were  required  to  provide  the  failure  pressures.  The 
failure  equations  for  lobar  buckling  and  general  instability  are  generally  from  reference  1, 
however  some  small  modifications  have  been  made.  A  complete  stress  calculation  is  included 
for  axisymmetric  yield. 

4.4.1  Axisymmetric  Yield 

The  majority  of  the  modifications  to  the  original  code  were  made  to  allow  calculation  of 
a  failure  pressure  for  this  mode.  The  beginning  code  provided  a  required  diving  depth  and 
therefore  a  pressure  to  withstand.  It  then  performed  a  von  Mises  stress  analysis  on  the  structure, 
resulting  with  a  shell  stress.  This  shell  stress  (after  modification  by  a  safety  factor)  was  then 
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compared  to  the  yield  strength  of  the  material  for  acceptability.  The  code  was  modified  to  iterate 
on  the  pressure,  with  the  failure  pressure  designated  as  the  pressure  at  which  the  calculated  stress 
was  equal  to  the  yield  stress  (without  any  modifications  by  safety  factors).  There  was  no  distinct 
equation  for  the  failure  pressure.  The  following  is  the  code  used  to  determine  osy,  which  is  the 
shell  stress. 

bv^p 

B:= - - — 

^  V  Area  ratio 


1 


cosh  (9)  -  003(0) 
sinh(9)  +  sin(0) 


Slenderness  Parameter 


Deflection  coefficient 


p:= 


2-N 

A  +  t^'tp 


Frame  flexibility  parameter 


Frame  deflection  parameter 


Bending  effect  (mem) 


He:=-2- 


3 

2 

1  -V 


sinh(0)  -  sin(9) 
sinh(9)  +  sin(0) 


The  following  begins  the  von  Mises  stress  analysis. 


Bending  effect  (bend) 


Bending  efffect  near  frame 
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4.4.2  Asymmetric  Buckling 

The  equation  for  lobar  bucking  is  the  same  as  that  used  by  SNAME  and  ABS  and  the 
analytic  solution;  it  is  the  Widenburg  approximation: 
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For  general  instability,  the  Bryant  equation  is  used  with  the  effective  length  of  shell  given 
by  the  following  method: 

Leff  W 

Where 

Lc  :=  Lf  -  tw  is  defined  as  the  clear  length  between  the  mid-bay  and  the  stiffener 

2  2 

4  cosh(ni-0)  -  cos(n2*0) 

0  cosh(n]’0)sinh(n}*0)  cos(n2’0)*sin  (r>2e) 
nj  n2 

nj  :=  0.5-^  1  -  y  n2  :=  0.5-/TTy 


F]  is  the  first  transcendental  equation  determined  by  Pulos  and  Salerno.  The  original 
methodology  required  the  pressure  as  an  input  to  first  determine  the  0  function  (it  is  the  same  as 
for  axisymmetric  yield)  and  then  Fi  was  found.  To  provide  for  a  simpler  method,  a  series  of 
pressures  that  encompassed  the  predicted  range  of  failure  were  input  into  the  above  equations 
and  an  average  value  of  Fi  was  chosen  from  the  resultant  values.  The  Fi  values  changed  by 
approximately  +  0.05  for  each  test  cylinder.  The  effective  length  was  then  used  to  calculate  the 
moment  of  inertia  for  the  combined  shell  and  stiffener.  The  number  of  lobes  was  set  to  be 
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between  2  and  4;  the  failure  pressure  was  then  found  without  other  modifications  to  the  Bryant 
formula. 


(22) 


Where: 


Chapter  5:  Experimental  Results 


The  various  Classification  Societies’  design  rules  were  tested  against  data  collected  from 
experiments  conducted  by  the  U.  S.  Navy  in  support  of  submarine  design.  Each  society’s  design 
rules  were  used  to  determine  the  failure  pressure  and  failure  mode  of  each  of  the  test  cylinders. 
The  resulting  predictions  were  then  compared  to  the  experimental  results.  Of  primary  interest 
was  agreement  between  the  design  rule  prediction  and  experiment  on  the  mode  of  failure, 
followed  by  the  closeness  of  the  predicted  failure  pressure  to  the  actual  collapse  pressure. 

5.1  NAVSEA  Test  Cylinders 

The  test  data  was  provided  by  the  Naval  Sea  Systems  Command  submarine  structures 
unit  (SEA  05P2).  Four  test  cylinders  were  selected  that  covered  the  range  of  examined  failure 
modes.  The  cylinder  diameter  to  thickness  ratios  (D/t)  fell  between  1 12  to  198,  modeling  typical 
submarine  D/t  ratios.  Two  of  the  cylinders  had  internal  stiffeners  while  the  other  two  cylinders 
had  external  stiffeners.  All  four  test  cylinders  had  built-up  end  stiffeners  with  a  combination  of 
slightly  different  spacing  and  /  or  larger  stiffener  dimensions  than  the  uniform  field  stiffeners. 
These  end  stiffeners  were  designed  to  prevent  shell  yielding  in  the  end  bays  due  to  increased 
stress  levels  associated  with  the  boundary  conditions.  It  was  estimated  that  without  the  end 
stiffeners  a  4-5%  reduction  in  axisymmetric  yielding  pressure  could  occur.  None  of  the  end 
stiffeners  met  the  classification  societies’  requirement  for  a  “deep”  stiffener.  None  of  the  design 
rule  codes  allowed  for  these  variable  stiffeners,  therefore  the  non-uniformities  were  disregarded 
and  the  end  stiffeners  treated  as  field  stiffeners.  The  four  cylinders  are  described  below. 
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5.1.1  Cylinder  l.d 


Cylinder  l.d  was  a  machined  cylindrical  shell  with  rectangular  external  ring  stiffeners. 
The  material  was  high  strength  steel  with  a  yield  strength  of  80,000  psi.  Figure  5  shows  the 
schematic  of  the  cylinder.  The  boundary  conditions  consisted  of  one  end  being  fully  fixed;  the 
other  end  had  all  freedoms  fixed  except  for  axial  displacement  (these  conditions  conflict  with  the 
design  rules  assumption  of  completely  clamped  ends).  External  hydrostatic  pressure  was  applied 
including  axial  line  load  to  simulate  load  on  the  end  plate.  The  experiment  tested  the  ability  of 
the  design  rules  to  predict  elastic  shell  bucking  (Lobar  buckling).  The  experimentally 
determined  collapse  pressure  was  633  psi  with  failure  by  asymmetric  (Lobar)  buckling. 
Appendix  A  has  the  analysis  of  predicted  failure  for  this  cylinder. 

Figure  5:  Test  Cylinder  l.d  Schematic 


44 


5.1.2  Cylinder  l.f 


Cylinder  1  .f  was  a  cylindrical  shell  with  internal  tee  stiffeners  of  welded  construction. 

The  material  was  high  strength  steel  with  a  yield  strength  of  98,500  psi.  The  boundary 
conditions  consisted  of  4.0  inch  steel  plates  attached  with  full  fixity  to  the  end  of  the  adaptor  ring 
on  the  model.  External  hydrostatic  pressure  was  applied.  This  test  cylinder  was  used  to  predict 
failure  by  elastic  general  instability.  There  was  no  experimental  elastic  collapse  pressure; 
therefore  the  critical  pressure  was  calculated  by  two  separate,  reliable  analysis  programs  with  the 
results  being  4858  psi  (with  3  waves)  and  4953  psi  (with  3  waves).  The  test  cylinder  actually 
failed  by  inelastic  general  instability  at  a  pressure  of  2200  psi.  Figure  6  shows  the  cylinder 
dimensions  and  Appendix  B  has  the  analysis  of  predicted  failure  for  this  cylinder. 


Figure  6:  Test  Cylinder  l.f  Schematic 


5.1.3  Cylinder  2.a 

Cylinder  2.a  was  a  machined  cylindrical  shell  with  external  tee  stiffeners.  The  material 
was  high  strength  steel  with  a  yield  strength  65,500  psi.  Figure  7  shows  the  schematic  of  the 
cylinder.  The  boundary  conditions  consisted  of  end  closures  made  of  3.0  inch  steel  plates 
attached  to  the  idealized  adaptor  ring  with  full  fixity.  External  uniform  hydrostatic  pressure  was 
applied  to  the  model.  This  cylinder  was  used  by  the  Navy  to  predict  end  bay  failure  (shell 
collapse  influenced  by  end  bay  design).  This  is  a  specific  example  of  axisymmetric  buckling  and 
was  used  as  the  axisymmetric  model  for  the  classification  society  rules.  The  experimental 
collapse  pressure  was  found  to  be  921  psi  by  axisymmetric  collapse  in  the  second  bay  from  the 
adaptor  ring.  Appendix  C  has  the  analysis  of  predicted  failure  for  this  cylinder. 

Figure  7:  Test  Cylinder  2.a  Schematic 
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5.1.4  Cylinder  2.c 


Cylinder  2.c  was  a  fabricated  cylinder  with  internal  ring  stiffeners.  The  base  material 
was  high  strength  steel  with  a  yield  strength  of  157,000  psi.  Figure  8  shows  the  schematic  of  the 
cylinder.  The  shell  was  cold  rolled  and  fabricated  with  a  deliberate  out-of-roundness 
imperfection.  The  frames  were  built-up.  The  frame  web  material  was  base  metal,  and  the  frame 
flanges  were  cold  rolled.  The  boundary  conditions  consisted  of  one  end  being  fully  fixed  with 
the  other  end  having  all  freedoms  except  axial  displacement.  External  uniform  hydrostatic 
pressure  with  an  axial  end  load  to  simulate  end  plate  loading  was  applied.  This  test  cylinder  was 
used  by  NAVSEA  to  predict  the  inelastic  general  instability  failure  mode  and  to  model  out-of- 
roundness  imperfections.  In  the  current  comparison  the  out-of-roundness  was  disregarded.  The 
collapse  pressure  was  experimentally  found  to  be  3640  psi  in  2  circumferential  waves  in  an 
inelastic  general  instability  mode.  Appendix  D  has  the  analysis  of  predicted  failure  for  this 
cylinder. 

Figure  8:  Test  Cylinder  2.c  Schematic 


5.2  Calculation  to  Experiment  Comparison 


Table  1  shows  the  comparison  of  the  classification  societies’  design  rules  to  the 
experimental  results.  The  table  displays  the  calculated  failure  pressure  and  failure  mode  for  each 
cylinder  and  the  percent  difference  from  the  experimental  failure  pressure  (if  the  experimental 
and  calculated  failure  modes  are  the  same). 


Table  1;  Comparison  of  Design  Rule  Calculations  to  Experimental  Results 


1  .d 

2. a 

2.C  1 

Pressure 

Mode 

Pressure 

Mode 

Pressure 

Mode 

NAVSEA 

EPERIMENT 

633 

H 

2200 

iGI 

921 

AX 

3640 

iGI 

Analytic 

Solution 

605 

-4 . 6% 

■ 

2141 

AX 

876 

-5.1% 

AX 

4080 

AX 

ABS 

605 

-4 . 6% 

■ 

2039 

AX 

844 

-8.8% 

AX 

4211 

AX 

PNA 

605 

“4 .6% 

■ 

1928 

AX 

815 

-12 .1% 

AX 

3864 

AX 

Germani s  cher 

Lloyd 

606 

“4.5% 

■ 

AX 

1030 

12 .4% 

AX 

4567 

AX 

13A  Professional 

Summer 

605 

”4 . 6% 

■ 

819 

-11 . 6% 

AX 

3712 

AX 

Key: 

L 

Asymmetric  (Lobar)  Buckling 

AX 

Axisymmetric  Yielding 

GI 

General  Instability 

e/i 

elastic  /  inelastic 

There  is  excellent  agreement  between  the  experimental  data  and  the  calculations  for 
cylinder  1  .d.  The  Lobar  buckling  failure  was  expected  as  the  slenderness  ratio  was  201 ,  several 
magnitudes  greater  than  the  breakpoint  of  1 . 1 4  between  asymmetric  and  axisymmetric  failure. 
The  agreement  between  the  design  rules  and  the  analytic  solution  calculations  was  expected  as 
they  all  used  the  same  formula  to  determine  the  lobar  buckling  pressure.  The  calculated  failure 
pressure  was  4%  below  the  experimental  pressure.  The  higher  experimental  critical  pressure  was 
attributable  to  the  test  cylinder  being  more  fully  clamped  than  theorized  in  the  design  rules. 
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Cylinder  2. a  also  generated  agreement  between  the  experimental  data  and  the  design  rule 
calculations.  The  predicted  failure  mode  of  axisymmetric  yield  was  confirmed  by  the 
experiment.  However  there  were  some  significant  differences  between  the  predicted  collapse 
pressure  and  the  experimental  failure  pressure.  All  of  the  design  rules,  except  for  the 
Germanischer  Lloyd,  predicted  failure  at  a  pressure  lower  than  that  experimentally  found.  The 
Germanischer  Lloyd  calculation  over  estimated  the  failure  pressure  and  was  the  furthest  from  the 
experimental  pressure  with  a  pressure  1 1 .8%  over  the  experimental  pressure.  This  failure  mode 
is  the  most  dependent  upon  the  yield  strength  of  the  material  therefore  a  small  variance  between 
the  given  yield  strength  and  the  actual  yield  strength  of  the  test  specimen  may  have  contributed 
to  the  differing  pressures  predicted  by  the  design  rules.  All  of  the  design  rules  use  a  simplified 
version  of  the  methodology  presented  by  Pulos  and  Salerno  [7].  The  Germanischer  Lloyd  and 
the  modified  Professional  Summer  calculations,  perform  an  iterative  operation  to  find  the 
pressure,  very  close  to  the  analytic  methodology.  The  ABS  and  PNA  calculations  use  single 
value  equations  substituted  for  the  transcendental  functions  of  the  analytic  solution. 

Cylinder  1  .f  failed  at  2200  psi  experimentally  in  an  inelastic  general  instability  mode. 

The  experiment  was  done  to  test  the  ability  of  NAVSEA’s  computer  codes  to  predict  elastic 
general  instability.  A  predicted  failure  pressure  of  4858  psi  was  determined  for  the  elastic 
general  instability  mode.  The  design  rule  codes  estimated  the  general  instability  failure  pressure 
very  well.  Table  2  compares  the  design  rules  and  analytic  solution  general  instability  pressures 
to  the  experimental  failure  pressure.  The  agreement  between  the  design  rules  was  expected  as 
they  all  use  the  Bryant  equation  (equation  10)  to  determine  the  failure  pressure.  The  only 
differences  come  from  the  variations  in  the  effective  length  of  the  shell  for  the  combined  shell 
and  stiffener  calculations  and  small  variations  in  the  radius  used  in  the  equation. 
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Of  further  interest  to  cylinder  1  .f  was  that  the  design  rules  predicted  failure  by 
axisymmetric  yielding  at  an  average  of  2206  psi..  As  the  cylinder  actually  failed  by  inelastic 
general  instability  the  design  code  predictions  indicate  that  the  two  modes  of  failure  are  very 
close  together.  This  cylinder  was  very  close  to  failure  in  multiple  modes  at  approximately  the 
same  pressure  (i.e.  an  identical  cylinder  of  the  same  dimensions  and  material  may  have  failed  by 
axisymmetric  yielding  vice  the  general  instability  depending  on  the  eccentricity  of  the  cylinder 
and  other  small  defects).  This  multiple  failure  mode  condition  must  be  guarded  against  in  real 
designs,  usually  by  applying  different  safety  factors  to  the  various  modes.  [1] 

Table  2:  Cylinderl.f  Elastic  General  Instability  Failure  Pressures 


DESIGN  RULE 

FAILURE  PRESSURE 
(PSI) 

%  FROM 
EXPERIMENT 

LOBES 

Experiment 

4858 

— 

3 

Analytic  Solution 

4496 

-7.5 

3 

ABS 

4496 

-7.5 

3 

PNA 

4496 

-7.5 

3 

Germanischer  Lloyd 

4651 

-4.3 

3 

13A  PS 

4460 

-8.2 

3 

Cylinder  2.c  failed  at  3640  psi  in  an  inelastic  general  instability  mode  with  2 
circumferential  waves.  However  the  design  rule  calculations  all  predicted  failure  by 
axisymmetric  yielding  at  an  average  pressure  of  4086  psi.  For  this  cylinder  the  design  rule  codes 
were  not  close  in  predicting  the  elastic  general  instability  failure  pressure.  This  large 
overestimation  of  the  failure  pressure  can  be  attributed  to  the  assumption  of  perfect  circularity  in 
the  design  codes,  whereas  the  experimental  model  had  a  two-wave  sinusoidal  imperfection  of 
maximum  height  of  +  0.105  inches.  This  deliberate  out-of-roundness  would  significantly  reduce 
the  resistance  to  buckling.  Table  3  compares  the  design  rules  /  analytic  solution  to  the 
experimental  failure  pressure. 
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From  Table  3,  the  axisymmetric  failure  pressures  were  higher  than  the  actual  failure 
pressure  but  all  were  within  16%  of  the  experimental  result.  This  closeness  between  the  failure 
modes  resembles  the  results  for  cylinder  l.f  Further  study  may  be  warranted  to  explore 
connections  between  axisymmetric  yield  and  inelastic  general  instability. 


Table  3:  Cylinder  2.e  Elastic  General  Instability  Failure  Pressures 


DESIGN  RULE 

FAILURE  PRESSURE 
(PSD 

%  FROM 

EXPERIMENT 

LOBES 

Experiment 

3640 

— 

2 

Analytic  Solution 

8642 

137 . 4 

2 

ABS 

8642 

137.4 

2 

PNA 

8642 

137.4 

2 

Germanischer  Lloyd 

9702 

166.5 

2 

13A  PS 

8536 

134.5 

2 
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Chapter  6:  Conclusions 

The  classification  society  design  rules  studied  are  important  tools  for  engineers  and  naval 
architects  designing  and  studying  cylindrical  structures  subjected  to  external  hydrostatic 
pressure.  The  engineer  must  have  confidence  that  the  design  code  used  will  provide  acceptable 
(safe)  calculations  for  his  or  her  structure.  This  confidence  can  be  assured  by  comparison  of  the 
calculated  failure  pressure  and  mode  to  that  found  from  experiments.  This  thesis  attempted  to 
provide  that  comparison  for  several  of  the  most  used  design  rules  along  with  a  comparison  of  the 
analytical  solution  upon  which  most  the  design  rules  are  at  least  partly  based. 

6.1  Comparative  Analysis  Review 

As  discussed  in  Chapter  5,  the  design  rules  had  mixed  results  in  correctly  predicting  the 
failure  mode  and  failure  pressure  of  the  test  cylinders.  In  general  all  of  the  design  rules  and  the 
analytic  solution  were  in  agreement  for  the  specific  geometries.  For  the  test  cases  of  failure  by 
axisymmetric  yielding  and  lobar  buckling,  the  calculated  pressures  were  accurate  when 
compared  against  the  experimental  results.  However  for  the  two  cylinders  that  experimentally 
failed  by  general  instability,  all  of  the  design  codes  predicted  failure  by  axisymmetric  yielding 
vice  the  general  instability.  The  design  rules  only  account  for  elastic  general  instability  which 
will  occur  at  a  higher  pressure  than  the  inelastic  general  instability.  For  cylinder  l.f  the  design 
codes  calculation  of  the  elastic  general  instability  mode  fmlure  pressure  was  very  accurate,  but 
for  cylinder  2.c  the  calculated  failure  pressures  were  all  greater  than  100%  over  the  experimental 
pressure.  The  majority  of  the  pressure  differentials  for  2.c  can  be  attributed  to  the  built-in  non¬ 
circularity  in  the  test  cylinder. 
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6.2  Agreements  and  Differences 


In  this  thesis  there  were  thirty  failure  pressures  calculated  (five  failures  for  each  of  four 
test  cylinders  and  five  general  instability  pressures  each  for  the  two  cylinders  that  failed  by 
general  instability)  with  design  rule  /  analytic  solution  agreement  on  33%  of  the  pressures.  The 
other  twenty  pressure  calculations  varied  from  over  30%  to  less  than  1%  different. 

The  pressure  calculations  that  were  in  agreement  were  generated  from  the  use  of  the  same 
equation  and  same  dimensions  for  most  of  the  asymmetric  buckling  and  elastic  general 
instability  calculations.  For  the  asymmetric  buckling  predictions,  ABS,  PNA,  and  the  MIT  1 3A 
Professional  Summer  calculations  use  the  same  equation  as  the  analytic  solution.  This  equation 
was  independent  of  the  number  of  circumferential  lobes  in  the  failed  part  of  the  shell.  For  the 
general  instability  cases,  the  ABS  and  PNA  calculations  agreed  with  the  analytic  solution.  These 
three  codes  used  the  same  effective  length  (Lc  =  frame  spacing)  and  the  assumption  that  the 
radius  of  the  combined  shell  and  stiffener  was  the  mid-plane  radius  of  the  shell. 

The  difference  between  calculations  can  be  accounted  for  individually:  The 
Germanischer  Lloyd  code  for  asymmetric  buckling  used  an  iterative  process  that  was  dependent 
of  the  number  of  lobes  in  the  failed  cylinder;  however  the  failure  pressure  was  still  within  a  tenth 
of  a  percent  of  the  other  calculations.  The  predictions  for  elastic  general  instability  were  also 
generally  close.  The  Germanischer  Lloyd  and  MIT  13A  Professional  Summer  codes  used 
different  assumptions  for  Lc  than  the  other  codes,  which  gave  predicted  failure  pressures  slightly 
different  than  the  base  general  instability  equation  (equation  1 0).  The  largest  variations  between 
the  codes  were  generated  for  the  axisymmetric  yield  mode.  All  of  the  codes  used  simplified 
variations  the  analytic  solution  developed  by  Pulos  and  Salerno.  These  variations  in  the  shell 
yielding  pressure  come  about  due  to  the  different  simplifications  made  in  the  codes. 
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6.3  Applications  of  the  Models 

The  various  design  rules  studied  are  promulgated  to  ensure  safe  design  of  cylinders  for 
use  under  external  pressure  conditions.  This  safe  operation  requires  an  almost  absolute  certainty 
that  the  design  will  not  fail  under  the  worst  anticipated  condition  (many  theories  and  practices 
exist  on  risk  based  design  and  the  use  of  safety  factors).  The  comparison  of  the  design  rules  to 
an  analytic  solution,  analysis  tools  (MIT  13A  Professional  Summer  and  SNAME),  and  to 
experimental  results  allows  a  designer  to  have  a  good  idea  of  how  a  particular  cylinder  design 
would  be  evaluated  by  each  entity.  This  comparison  ability  would  be  useful  in  judging  the  initial 
feasibility  of  a  design  and  also  would  be  useful  in  applications  where  a  design  would  be  subject 
to  more  than  one  classification  society. 

The  comparisons  in  this  thesis  should  in  no  way  be  used  as  a  detailed  design  tool  for  the 
subject  cylinders.  After  an  initial  design  is  compared  and  judged  to  be  adequate,  much  more 
rigorous  analyses  must  be  used  to  ensure  a  safe  design.  These  advanced  analyses  should  include 
finite  element  methodology  and  other  tools  that  can  look  at  local  stresses  instead  of  generalized 
stresses  in  the  shell.  These  higher  order  analysis  tools  can  account  for  material  differences, 
geometric  eccentricities  (out-of-roundness  and  other  along  the  shell),  varied  spacing  and  sizes  of 
stiffeners,  and  actual  construction  factors  such  as  heat  effected  zones  around  welds,  and 
bulkhead  effects.  Another  important  area  that  is  addressed  by  other  parts  of  the  design  rules  but 
not  studied  here  is  that  of  shell  penetrations.  These  discontinuities  in  the  shell  are  very 
susceptible  to  stress  concentration  and  must  be  reinforced  to  prevent  failure  at  pressures  lower 
than  that  predicted  for  a  continuous  shell. 
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6.4  Further  Areas  of  Study 

There  are  several  areas  that  require  further  research  to  completely  understand  the  result  of 
the  above  analysis.  The  area  most  evident  in  need  of  more  study  is  the  failure  of  the  cylinders  by 
inelastic  general  instability.  While  some  work  has  been  accomplished  on  this  phenomenon,  most 
of  the  work  is  either  classified  by  government  entities  or  is  empirical  data.  As  part  of  the 
research  into  the  inelastic  general  instability  should  be  a  detailed  analysis  of  interactions  / 
relationships  with  axisymmetric  yield.  As  found  from  the  experimental  cylinders  that  failed  by 
inelastic  general  instability,  the  design  codes  and  analytic  solution  predicted  failure  by 
axisymmetric  yield  of  the  shell  at  pressures  close  to  the  actual  failure  pressures.  These 
comparisons  suggest  that  at  least  the  two  modes  are  very  close  together  and  may  have  some 
interaction. 

There  are  many  other  classification  societies  that  produce  design  rules  for  stiffened 
cylinders  and  other  geometries.  These  societies  include  the  American  Petroleum  Institute  (API), 
NORSOK  and  DetNorske  Veritas  (Norway),  Lloyd’s  Register  (United  Kingdom),  Registro 
Italiano  Navale  Group  (RINA)  (Italy)  and  several  others.  These  additional  design  rules  could  be 
compared  against  the  existing  test  cylinders  as  well  as  all  of  the  design  rules  should  be  compared 
against  more  experimental  data. 

Finally,  the  base  methodology  can  be  expanded  to  include  other  geometries  such  as 
spheres,  hemi-spheres,  cones  and  toroids.  The  classification  societies  address  most  of  these 
forms  which  can  be  found  as  the  end  closures  for  most  submersibles  and  other  pressure  vessels. 
These  more  complex  geometries  present  more  challenging  analytic  solutions  and  are  in  general 
harder  to  manufacture.  Experimental  data  is  scarcer;  therefore  finite  element  models  may  be 
needed  to  generate  comparisons. 
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Appendix  A:  Codes  for  Test  Cylinder  l.d 
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AMERICAN  BUREAU  OF  SHIPPING 


Rules  for  Building  and  Classing  Underwater  Vehicles, 


Systems,  and  Hyperbaric  Facilities 

Definitions 

ksi  :=  lOOOiJsi 

E  :=  3000(ksi 

Modulus  of  Elasticity 

V  ;=  0.3 

Poisson's  Ratio 

Gy  :=  8000Q)si 

Yield  strength 

Shell  Parameters 

Lg  :=  4.266n  Distance  between  stiffeners 

=  22.488in  Distance  between  bulkheads 

R  =  8.007in  Mean  radius  of  shell 
t  =  0.081  in  Thickness  of  shell 

R„  ;=  R  +  -  Outer  radius  of  shell 

o  2 

Dq  :=  2  Rjj  Outer  Diameter  of  shell 

Ring  Stiffeners 

t^;=0.l38n  thickness  of  web  of  ring  stiffener 

depth  :=  0.5'^n  height  of  ring  stiffener 

b  :=  faying  width  of  stiffener  (from  P&S  for  I  beam  stiffener) 

bf  :=  o.Oin  breadth  of  ring  stiffener 

b2  :=  bf  -  breadth  of  ring  stiffener  minus  the  web  thickness 
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d  :=  O.Oin 


ring  stiffener  flange  thickness 


L|,  :=  Lg  - Unsupported  spacing  between  stiffeners 
L:= 


L  =  4.266in 


C]  •=  -S- 


depth  +  b2d 
^  tv^-depth  +  bjd 


first  centroidal  height  of  ring  stiffener 


C|  ==  0.285in 

C2  :=  depth  -  cj  second  centroidal  height  of  ring  stiffener 


C2  =  0.285in 

h  :=  C]  -  d  distance  from  centroid  of  ring  stiffener  to  nearest  edge  of  flange 
Rg  ;=  R  +  .5  t  +  C2  radius  to  centroid  of  ring  stiffener 
Rg  =  8.332in 

Ag  :=  depth  +  b2  d)  cross-sectional  area  of  ring  stiffener 


Ag  =  0.079in^ 


:=  f-Y^bf-C]^  -  b2h^  +  Vc2^^  moment  of  inertia  of  ring  stiffener  about  its  centroidal  axis 


I^=  2.13x  10 


Rf  ;=  R  +  .5  t  +  depth  Radius  to  tip  of  the  stiffener 
R^=  8.61 7in 

Combined  Plate  and  Ring  Stiffener 


:=  mini 


i.s-VrI 

V  0-75Lg  J 


effective  length  of  shell  plate 
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area  of  effective  plate 


Lg  =  1.208in 

Ap:=Lt 

Asp  :=  Ap  +  As  area  of  plate  and  ring  stiffener 

Hj,  :=  depth  +  t  height  of  combined  plate  and  ring  stiffener 

Hj.  =  0.054a 

Bj  :=  L  - plate  length  minus  the  web  thickness 

111 

B,t^+  b2d(2-H^,-d) 

‘^ic  •=  - n — I - 77 -  neutral  axis  of  combined  plate  and  ring  stiffener 

+  Bjt  +  002 

from  outer  fiber  of  plate  (Rq) 

C|^  =  O.lOlin 

:=  R  -  .5 1  +  C|^  radius  to  centroidal  axis  of  combined  ring  stiffener  and  shell 
R^  =  8.067in 

•e  •=  -  (bf  -  t  J  (Hj,  -  cj^  -  d)^^ 

moment  of  inertia  of  combined  plate  and  shell 

Ig  =  O.OOOin'* 


General  Equations 


M  =  5.297 
1 
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e  =  6.809 


Q  =  3.405 

^  cosh  (2- Q)  -  cos(2  Q) 

N  := - 

sinh(2-Q)  +  sin(2*Q) 

N  =  0.997 

^  ^  (sinh(Q)  cos(Q)  +  cosh(Q)  sin(Q)) 

sinh(2Q)  +  sin(2Q) 

G  = -0.081 


^  _  sinh(2Q)  -  sin(2  Q) 
sinh(2-Q)  +  sin(2Q) 

H  =  0.998 

Inter-Stiffener  Strength  (6.19.1) 

1)  Inter-stiffener  strength  equations 


This  equates  to  axisymmetric  bucklint 


Effective  area  of  plate  and  stiffener  (External  stiffeners) 


A  =  0.073in^ 


F;= 


A  +  + 


2NtL 

0 


F  =  -0.027 


yield  pressure  at  midbay  and  midplane  of  cylinder 
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Py  =  787.9psi 


This  corresponds  to  Lobar  bucklinr 


von  Mises  buckling  pressure  for  a  cylinder 


(.-T. 


—  -  0.4sf  — 
2-R  l2-R 


=  604.9psi 


Maximum  allowable  working  pressure  for  inter-stiffener  strength 


P  P 

p  :=  if  JI1<, 

2  P 

y 

p.  |_^  if  1<>S3 


5 

-  Py  if  -  >3 

6  y  Py 


—  =  0.768 


P«  -  302.4psi 


^aits  * 


Paits  =  242psi 
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2)  longitudinal  stress 

Limiting  pressure  corresponding  to  the  longitudinal  stress  at  stiffeners  reaching  yield.  No  direct 
correlation  to  major  failure  modes 


A-  1  -  - 

I  2 


A  +  t^t  + 


2-N-t-L 


f  12  r 

P,:= - -  1+  -  y-H 

'  R  2 

V  ]  -  V  / 


P|  =  732.7psi 


Pals  :=  Pi-.67  Maximum  allowable  working  pressure  for  longitudinal  stress 


Pgjs  =  490.88psi 


b  Overall  Buckling  Strength  (6.19.5) 


A2(n)  :=  n  -  1 


A,(n)  := 


!]-(n^  .  xf 


Et  E'[gA2(n) 

Pnl(")  :=  ^ - 

^  R^-L 
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Given 

n  >  2 

n2  :=  Minimiz^Pj^  j  ,n) 
n2  =  3.109 

n  must  be  an  integer 

n2int:=  round  (n2.0) 


"2int  -  ^ 

"2flrst  ~  "2int'^®®  '^'’2int”^ 

n2int  “  1  otherwise 

"2first  =  2 

"2first 
"2int 
"2int  +  ^ 

"2int  +  2 

"2prime,  =  ^ 


5.127X  10^ 
1.522X  10^ 
I.974X  10^ 
3.022X  10^ 


psi 


Pn  :=  min(p„2) 


66 


Pj^  =  1522.3psi 
Paol  ■=  Pn--5 
P30l  =  761.131psi 

Summary 

Axisymmetric  Buckling 
Py  =  787.9psi 
Lobar  Buckling 
Pjj^  =  604.9psi 
General  Instability 
P3  =  1522.3psi 

"2int  =  3 


Germanischer  Lloyd 


1988  Edition 

Chapter  2  -  Submersibles 

MPa  :=  lO^Pa  ksi  :=  lO^psi 

General  Definitions 

E  :=  3000Qcsi 


L|  :=  4.266n 

Stiffener  spacing 

-  22.48an 

length  of  cylinder  between  bulkheads  or  lines  of  support 

s  :=  .081in 

thickness  of  shell 

R  :=  S.OO^^n 

radius  to  centerline  of  shell 

V  :=  0.3 

Poisson's  Ratio 

k:=  8aksi 

minimum  yield  stress  of  material 

Stiffener  Dimensions 

:=  O.Oin 

flange  thickness 

d^:=0.0in 

width  of  flange  from  web  to  edge  of  flange 

b2:=  2-df 

d^  :=  0.5'^n 

Height  of  web 

:=  0.138n 

web  thickness 

b  ;= 

width  of  stiffener  ring  in  contact  with  shell 

Aj  :=  ef-(e^  +  2  df)  +  cross-sectional  area  of  stiffener  ring 

Aj  =0.079in^ 
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^^stifTF'=‘^w+  H 


Hstiff=0-57in 


2  2 
2  V^stiff  +  '^2'^f 

2  ^w'^stiff  + 


C|  =  0.285in 


^2  '-^stiff  “‘^1 


C2  =  0.285in 


’’2  •-  ~ 


h2  -  0.285in 


■]{K  +  2-df)-c/  -  (2-df).h2^  +  e^-C2^_ 


Centroidal  Moment  of  inertia  of  ring  stiffener 


Ij  =2.13x  10 


e  :=  C2  +  .5  s  distance  from  stiffener  centroid  to  center  of  shell 


e  =  0.325in 


Effective  Stiffener  and  Shell 


^etest  ■  ^  ^  yjl  K-s 


Lg  :=  Lgjgjj  if  Lgjggj  <  L]  Effective  length  of  shell  (eqn  45a, b) 

Lj  otherwise 


=  1.277in 


I  *"  ^1  ■* - 

®  A,  ’  12 


Moment  of  Inertia  of  combined  plate  and  shell 


Ig  =  6.92x  10 
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^stiffc  *  Total  height  of  stiffener  and  plate 


Hstiffc  =  0-65  lit! 


d,  :=s 


bs  1 .  Lg 


2  2 

2  ^w’^stiffc  +  ^2‘^f  “  ^) 


Cj^  =  0.47in 


^2c '  ^stiffc  ^Ic 


C2c  =  O.lSlin 


Rq  :=  R  -  .5  s  +  C2c  radius  of  stiffener  ring  centroid  including  effect  of 


R^  =  8.148in 


A  :=  Aj - Modified  area  of  stiffener  ring 


A  =  0.076in 


ei:=C2c  +  .5s 


Cj  =  0.222in 
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2.4  Asymmetric  Buckling  (Lobar  but  not  named  this) 


n  :=  2  initial  guess 

1  + 

12  R^  (l  -  v^) 

■  '  2  2 
n  —  1  +  .5’X| 

E-sP„](n) 

Pnl(")  := - - 

Given 

n  >  2 

Hg  Minimiz^Pj^j  ,n) 

==  10.612  n  must  be  an  integer 
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"2  :='^'n'm'^Pn]a>") 


n2  =  3.247 


n  must  be  an  integer 


"2int=3 

"2int  -  ^  ^ 

"2int 
"2int  +  ^ 

^2int  J 

Pnla^*^2prime^^ 

Pnla^^2primej  j 

Pnla^^2prim^^  ^ 

Pnla^^2prime^^ 

Pg  :=mir(pg2) 

Pg  =  1244psi 

2.6  Symmetric  Buckling  (axisymmetric) 

More  Definitions 
For  Elastic-Plastic  Region  (Eqns  are  in  the  program) 

Ej  :=  E  Secant  Modulus  =  Young's  Modulus  for  elastic  region 

Ej  :=  E  Tangent  Modulus  =  Young's  Modulus  for  elastic  region 

Vp  :=  V  Poissons  Ratio,  elastic-plastic;  =  Poisson's  ratio  in  elastic 
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a  =  0,063 - 

mm 


C5  :=  a-Lj  calculation  factor  for  symmetric  buckling 

C5  =  6.809 

2  s^  E  .  . 

Past  •= - — T critical  pressure,  elastic,  calculation  factor 

Past  =  3716.21  Ipsi 

C5  :=  a-Lj  calculation  factor  for  symmetric  buckling 

C5  =  6.809 

G:=-^ 

Past 

+  G 

'  2 

Cg  :=  Cj-Cg 
C9  :=  CjC-j 


cosh( 

Cg)' 

2  / 
-  cost 

C9) 

2 

cosh  I 

[cs] 

I'Sinh 

(Cg 

)  COS 

(C9 

)-sin| 

(c,) 

^7 


cosh^Cgjsin^Cp)  sinh^Cg^cos^C^) 


C 

7 

C 

6 

coshj 

(c*) 

l-sinhj 

(c*) 

1  cos| 

N 

l-sinl 

(^9) 

^7 
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cos^C^^sin^C^)  cosh(Cg)sinh^Cg) 


- 


2  cosh^Cg^  sinh^Cg)  cos  sin 


1  -  V 


cosh  ♦  sin^C^^  sinh  (Cg)  ■  cos 


F4:= 


^  2  cosh  (Cg^- sinh  (Cg)  cos  sin 


^lO’-’ 


1  -- 

2J 


s*L| 


+  —  + 


sLj  L| 


V  '■v 


C],  :- 


0.91 


1  -  V 


:=■ 


-pR 


'^(t)--^o  (^  “  Cio'p2  +  ''■^io'Cirf'4) 


'^xi'-%i  2  “ ^10'^ir^4 


ax=-3.286x  10  Pa 


'  “^10p2“'''^10'^irp4) 


ct^  =  -6.651x  10  Pa 


Ko:=- 

^x 


.8k  =  64ksi 


Determination  factor  (if>.8k  then  must  use  Et,  I 


•S,  etc) 
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Oj  =  83.538ksi 


a  =  96.5ksi 


E,  =  3.215X  10 


a  :=  k’  .8  +  .2  tanh  5- — £•  -  4 


a  =  79456psi 


k  /  f  „ 

£  := -•  0.8+ 0.2  atanh  5— 
E  I  Ik 


£  =  3.215X  10 


k  f  (  E 

Ej  ;=  — •  0.8+  0.2-  tanh  5 — e 

e  _  V  V  k 


Ej  =  24548.3  Iksi 


E,  ;=  E-  1  -  tanh  5 — £  -  4 
V  V  k 


E,  =  1019.96a<si 


>  M  ^  *^S 

V  ;= - - V  — 

E  2  l2  J  E 


Vp  =  0.336 


Paa  ■= 


trial2:=  +  0 

triaI2=  I 


iterate  on  m  until  trial  is  <  trial2 


Summary 

Axisymmetric  Buckling 
=  923.8psi 

Lobar  Buckling 
Pj,j  =  605.9psi 
General  Instability 
Pg  =  1244psi 
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Society  of  Naval  Architects  and  Marine  Engineers 


Submersible  Vehicle  Systems  Design  (1990) 
ksi  :=  100C[3si 

Input  Section 


Lf;=  4.266n 

Length  between  frames 

L|j  :=  22.48an 

Length  between  bulkheads 

R  :=  8.00'Jin 

Radius  of  cylinder  to  centerline  of  shell 

t:=  0.08  lin 

shell  thickness 

1=  R  +  .5’t 

radius  to  outside  of  shell 

R0=8 

.047in 

E;=  30000ksi 

modulus  of  elasticity 

Gy  :=  80- ksi 

minimum  yield  stress  of  material 

H  :=  0.3 

Poisson's  ratio 

Ring  Stiffeners 

t^:=0.13an 

w 

thickness  of  web  of  ring  stiffener 

depth  0.57in 

total  height  of  ring  stiffener 

bf:=  O.Oin 

breadth  of  ring  stiffener 

d  :=  O.Oin 

ring  stiffener  flange  thickness 

b2  :=  bf  - breadth  of  ring  stiffener  minus  the  web  thickness 
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first  centroidal  height  of  ring  stiffener 


2  2 
t^^,depth  +  b2d 

depth  +  b2d 
Cj  =  0.285in 

C2  •=  -  C]  second  centroidal  height  of  ring  stiffener 

C2  =  0.285in 

^  distance  from  centroid  of  ring  stiffener  to  nearest  edge  of  flange 

h  =  0.285in 

Rj  :=  R  +  .5  t  +  C2  radius  to  centroid  of  ring  stiffener 

Rr=  8.332in 

^r'“  +  b2  d)  cross-sectional  area  of  ring  stiffener 

=  0.079in^ 

~  *w‘=2^)  moment  of  inertia  of  ring  stiffener  about  its  centroidal  axis 

If  =  2.13 X  IO~^n'* 

Combined  Plate  and  Ring  Stiffener 

Ap  :=  Lf-t  area  of  effective  plate 

'^sp  "^p  +  plate  and  ring  stiffener 

:=  depth  +  t  height  of  combined  plate  and  ring  stiffener 

•=  k  “  ‘w  P'ate  length  minus  the  web  thickness 
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neutral  axis  of  combined  plate  and  ring  stiffener 


VH,^  +  Bj-t^  +  b2d(2H,-d) 

VHc  +  Bi-t  +  db2 

from  outer  fiber  of  plate  (RO) 

Cjj,  =  O.lOlin 

R^.  :=  R  -  .5-t  +  cjj.  radius  to  centroidal  axis  of  combined  ring  stiffener  and  shell 
Rj,  =  8.067in 

'e  ‘ ‘  - 'Ic  -  jf] 

moment  of  inertia  of  combined  plate  and  shell 

Ig  =  O.OOPin'* 

Axisymmetric  yielding 


e  =  6.809 

cosh  (9)  -  cos(0) 
sinh(9)  +  sin(9) 

N  =  0.997 


(5  =  1.126 
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Py  =  781.382psi 


Lobar  Buckling 


pjj  =  604.9psi 


General  Instability 

tc-R 
m  := - 

Lb 

n  :=  1 


Given 

n  >  1 
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n2  :=  Minimiz^Pj,^j,n) 

n2  =  3. 1 09  n  must  be  an  integer 

"2int  ~  ^ 

"2first-=  "2inf'®®  "2int=  '  "2flrst  =  2 

”2int  ~  ^  otherwise 

"2first 
"2int 
"2int  +  ^ 

"2int  +  2 

Pcrl(^"2prim^^ 

Pcrl^*’2primej^ 

Pcrl^”2prim^^ 

Pcrl^''2prim^^ 

Pcr:=min(pcr2) 

Pj.^  =  1 522.3psi 

Summary 

Axisymmetric  Buckling 
Py  =  781.4psi 
Lobar  Buckling 
P|j  =  604.9psi 
General  Instability 
Pj.j.  =  1522.3psi 
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MIT  13A  Professional  Summer  Submarine  Design  Trends 

PROGRAM  TO  COMPUTE  SUITABILITY  OF  SUBMARINE  DESIGN  PARAMETERS. 
Ref:  "Hull  Material  Trade  Off  Study",  D  Fox,  Jan  94 

ksi  H  1  OOQjsi 

Global  Variable  Inputs: 

esO.Oin  Eccentricity 

Material: 


Gy  =  SOOOQ^si 

yield  stress  of  material 

E  =  3000(ksi 

young’s  modululs 

V  =  0.3 

poisson’s  ratio 

Geometry: 

R  =  8.00'^n 

shell  radius 

Dh2R 

shell  diameter 

Lf  =  4.266n 

frame  spacing 

Ls  =  22.48an 

bulkhead  spacing 

tp  =  0.08  lin 

shell  thickness 

tf  =  O.Oin 

flange  tickness 

Wf  =  O.Oin 

flange  width 

=  0.13Sn 

web  thickness 

h^y  H  O.STin 

web  height 
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Compute  areas: 


R,  =  R  +  —  Frame  Radius 
2  . 

Rj  =  8.047in 

Af  ;=  tf  wf  Frame  flange,  web  area 

'^w  tw’hw 

A  ;=  Af  +  Aw  Frame  Area  =  Flange  +  Web 
A  =  0.079in^ 

PARTI  SHELL  YIELDING 

Von  Sanden  and  Gunther  (1952) 

PNA  Section  8.4 

Uv’tp 

B:- - — 

^  ^w-tp  Atrea  ratio 


B  =  0.124 


0  :=  Lf 


b1i-v^) 

(R-tp)^  , 


Slenderness  Parameter 


0  =  6.809 


cosh(0)  -  cos(0) 

sinh(0)  +  sin(0)  Deflection  coefficient 


p:=- 


2-N 

r  1  1 

0.25  r 

■  J 

A  +  t^tp 

1 

1 

< 

1 _ 

P  =  1.126 


Frame  flexibility  parameter 
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-B 


Frame  deflection  parameter 


r  =  0.341 


■J®  ufei  re 

smh  —  cos  —  +  cosh  —  -sin  — 

Uy  U;  UJ  I2 


sinh(0)  +  sin(0) 


Bending  effect  (mem) 


Hm  =0.08! 


He:=-2 


•  kI  ®  1  ( ®  1  ufe^  .  fe 

J_1  V2j  [2J  (2J  I2 


sinh(0)  +  sin(0) 


Bending  effect  (bend) 


He  =  0.085 


P  crl  =  f  y  <-  O  y 
P|,  1  psi 


sinh(0)  -  sin(0) 
sinh(0)  +  sin(0) 


Pj  ^  1000  psi 

P2  ^  2000  psi 

delta  <-  5  psi 
limit  1  psi 
conv<-  1  psi 
j^O 

while  j  <  20 
I  for  i €  0  .2 


-r .  rt. 

- [1  +  r  tHfvi  +  V  Hei] 

-P.  R 

V — +  r  iHM“  V  He'] 

I  -J 


-P.R 

^  xxso  ^ - +  F  H  E' 


Bending  efffect  near  frame 


-P.R 

®  xxsi  ^ - '0  5  -  F  H  E» 

-^i^\  [  f  3  ■ 

c«fo< - 1  -r-  1  +  V  -  K 

L  L  l.-v^J  . 

\  (  3  1 

OiWfi< - 1-r  i-v  -  K 

‘P  .  ll-v'J 
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Perl  ~ 


801.682 

14 


Peril  •"  Perlji'  ^P*' 

Peril  =  801.7psi 

PART  2  LOBAR  BUCKLING 

Windenberg  Approx  of  Von  Mises  (1933) 
Assumes  n  lobes  =  Pi*D/L 
Collapse  pressure: 


PcLB  =  604.9psi 

PART  3  GENERAL  INSTABILITY 

Corrected  Bryant  Formula  (1954)  for  better  model  test  correlation 
Pressure  loading  is: 

P;=  pgDfSFgi 
P  =  6697.99psi 

Compute  effective  frame  spacing: 


Y=  1.802 

Compute  clear  length: 
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Lq  Lf 

n\  :=  O.S-yJ  1  -  Y 

ni  0.448i 
Web  thickness: 
n2  :=  O.S>/ 1  +  y 

n2  =  0.837 

2  2 

4  cosh(ni*0)  -cos  ("2-6) 

0  cosh(n|-0)sinh(n|-9)  cos(n2-0)-sin(n2-0) 
ni  nj 

Fi  =  0.23 
must  be  less 
than  1 .00 

FI  is  almost  a  linear  decreasing  function  for  pressures  from  1  to  2000  psi  with  an  average  value 
of  0.27.  This  will  be  used  in  the  following  analysis  as  the  pressure  is  the  unknown  and  therefore 
the  above  equations  cannot  be  directly  used. 

Fi  :=  0.27 

Lgff  :=  Lc-Fi  +  tvv  Effective  shell  plate  length: 

Lgff  =  1 .253in 

Theoretical  critical  lobe  number  values  are: 


i:=0..2 


Circumferential  Lobes 


Aeff  :=  Lgff  tp  Effective  plate  area 
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Agff  =  O.lOlin^ 


m:=7t —  Longitudinal  Lobes 


m=  1.119 


hw  +  tf]  fh,v  +  tp^ 

.Af - —  -Acff 

^  /  V  2  y 

Aeff  +  A„,  +  Af 


Frame-plate  neutral  axis  (ref  web  centre+  toward 


flange); 


Yna  =  -0.183m 


Uses  Parallel  Axis  Theorm:  Icor  =  I  +  Ad2 


Moments  of  inertia  for  plate, flange, web: 


1  _  I  A  *p  +  '’w 

^pcor  --  Ip  ^eff  I  +  yna 
V  2  / 


I\vcor  •“  I\v  ■*"  ^w’(yna) 


1  _  T  .  A  If 

Ifcor  •-  If  +  Af  ~  yj^a 

V  2 


Total: 
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•“  %cor  Kvcor  +  ^fcor 


.  4 

Igff  =  0.007in 

The  critical  Elastic  General  Instability  pressure  is: 


Min  Pressure: 

f  3 

5.035X  10 

PcGl=  1,278x  10^ 

Vl.516x  10^ 

PcGI  :=  miri(PcGl) 

PcGI  =  I278.1psi 

Summary 

Axisymmetric  Buckling 
Peril  =  801.7psi 
Lobar  Buckling 
PcLB  =  604.9psi 
General  Instability 
PcGI  =  1278.1psi 
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Analytic  Solution 


Definitions 

ksi  :=  6.894757 loSa 

E:=  3000a<si 

Modulus  of  Elasticity 

R  :=  S.OO^n 

Radius  of  cylinder  to  centerline  of  shell 

D:=2  R 

Diameter  of  cylinder 

Gy  :=  80ksi 

Yield  strength 

L:=4.266n 

Length  of  supported  cylinder 

:=  22.48an 

Distance  between  bulkheads 

M  :=  3 

Poison's  ratio  for  Fe/Steel 

t  :=  0.08 lin 

Shell  thickness 

Ring  Stiffeners 

:=  O.JSSn  thickness  of  web  of  ring  stiffener 
H  :=  o.57in  height  of  ring  stiffener 

b  :=  faying  width  of  stiffener  (from  P&S  for  I  beam  stiffener) 

bf  ;=  o.Oin  breadth  of  ring  stiffener 

b2  :=  bf  - breadth  of  ring  stiffener  minus  the  web  thickness 

d  :=  O.Oin  ring  Stiffener  flange  thickness 

hjnash  :=  L  -  bf  distance  from  flange  edges 


92 


tw'H  +  bjd 


first  centroidal  height  of  ring  stiffener 


Cj  =  0.285in 


C2:-H-c, 


second  centroidal  height  of  ring  stiffener 


C2  =  0.285in 


h  :=  Cl  -  d 


distance  from  centroid  of  ring  stiffener  to  nearest  edge  of  flange 


h  =  0.285in 


Rj.  :=  R  +  .5  t  +  C2  radius  to  centroid  of  ring  stiffener 


R^=  8.332in 


:=  (t^  H  +  b2  d)  cross-sectional  area  of  ring  stiffener 


=  0.079in 


'X  'X  ^ 

Ij.  :=  -  l  l  bf  cj  -  b2  h  +  VC2  j  moment  of  inertia  of  ring  stiffener  about  its  centroidal  axis 


I^  =  2.13x  10 


yff-=y 

I  R^ 


effective  area  of  stiffener  eqn  [24a]  from  P&S 


Agff  =  0.076in 


a  := -  ratio  of  effective  frame  area  to  shell  area  eqn  [62]  P&S 

Lt 


a  =0.219 


p  :=  —  ratio  of  faying  width  to  frame  spacing  eqn  [62]  P&S 


P  =  0.032 
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Combined  Plate  and  Ring  Stiffener 

Ap  ;=  L  t  area  of  effective  plate 

Ajp  :=  Ap  +  A^  area  of  plate  and  ring  stiffener 

H(,  :=  H  +  t  height  of  combined  plate  and  ring  stiffener 

Hj,  =  0.054ft 

B,  :=  L  - plate  length  minus  the  web  thickness 
VH/+B,t^+b2d(2H^-d) 

C]c  :=  -5 -  neutral  axis  of  combined  plate  and  ring  stiffener 

+  B|  t  -f  d*b2 

from  outer  fiber  of  plate  (Rq) 

c  =  O.lOlin 

:=  R  +  .5  t  +  radius  to  centroidal  axis  of  combined  ring  stiffener  and  shell 

R^,  =  8.148in 

'e  -  Br(c,c  -  f  +  bf.(H,  -  -  (bf  -  g.(H,  -  c,,  -  d)^] 

moment  of  inertia  of  combined  plate  and  shell 
Ig  =  O.OI  in^ 

Buckling  of  unreinforced  shells 

von  Mises  bucking  pressure: 

n:=  1 

guess  at  number  of  waves  around  the  circumference 
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Axisymmetric  Buckling 


P  c  a  b  * 


y  0 

1  im  it  I  0  p  si 
t  est  <—  0  p  si 
c  0  n  V  <—  I  p  si 

j  ^  0 

wh  i  le  j  ;<  2  0 

I 


T|  ,  ♦-  -  yj]  -  Y 


I 


T|  2  <-  j  yj  1  +  y 

e  -1/3  (1  -  n^).— ^ 


cosher)  -  cos^il  2  ®)^ 


G  cosh^rij'Gysinh^r|j*0^  cos^Tl2G))'Sin^r|2'0) 

1  ^2 
cosh^T'i  J-G^  sin^ri  2-®)  sinh^ri  j-Gycos^Tl  2'®) 

_ _ 

cosh^n  |-0)sinh(r|  jO^  cos(r|  2  ^^)sin  ^t|  2-0^ 

n  1  n  2 


<  1  ^ 

^  2 

^  1 

J  1  -  Li"'  (’1 

J  O^  sin  h  (r|  |  0^ 

cos(r|  2-G)'sin  ^ri  2'^) 

n  1 

^  2 

c  0  sh  ^r| 

I  *0^ sin  2  '^ ) 

sinh^Tj  |•0^■cos^r^  2*^) 

,  ,  L± . 

n  2 

^  1 

J  1  - 

j  ‘0^ sin  h  ^ri  j  -0^ 

cos(ri  2  0)  sin  3 

^  I 

n  2 

0-f) 


«  +  p  +  (1  -  p)-F  J 

de  n  o  m  1  < —  A  -j  F  2  "*■  2*^ 


dc  n  o  m  2 


P  c2 


ii 


+  F4^{l-|a-f^^ 


*4-  dcnom  1  “denom2 


break  i  f  |  p  ^  2  "  *  |  <  I  im  it 

.est  t-  p  ,2 

j  <-  j  I 

P  c2 
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(  885.4^ 
Pcab  =  [  2  j 

Paxiy“  PcabQ'^P^' 


Paxiy 


885.4psi 


Asymmetric  Collapse  (Lobar  buckling) 

Windenberg  Approx  of  Von  Mises  (1933) 
Assumes  n  lobes  =  Pi*D/L 


PcLB=  604.9psi 


General  Instability  of  shells  and  rings 


Given 

n  >  1 

n2 Minimiz^p^Qj,n) 

n2  -  3.109  n  must  be  an  integer 
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n2int;=round(n2,0) 


"2int  -  3 


"2p: 


nme*“ 


"2int  -  ’ 
"2int 
"2int  +  ' 
"2int  +  2 


Pcgi2  ■■ 


) 

PcGl^*’2prime|j^ 

PcGl^’’2primej^ 

PcG^*’2prime^^ 

PcGl^'’2primej^ 


Pcgibryant  •“  '’^^Pcgi2) 


Pcgibryant  =  ••522x  10  psi 


Summary 
Lobar  Buckling 
PcLB  =  604.9psi 
Axisymmetric  Buckling 
Paxiy= 

General  Instability 
Pcgibryant  ~  1522.3psi 
"2int  ~  ^ 
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Appendix  B:  Codes  for  Test  Cylinder  l.f 
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AMERICAN  BUREAU  OF  SHIPPING 


Definitions 

ksi  :=  100C|3si 

E  :=  3000(ksi  Modulus  of  Elasticity 

V  :=  0.3  Poisson's  Ratio 

Oy  :=  9850q3si  Yield  strength 

Shell  Parameters 

Lj  :=  2.665n  Distance  between  stiffeners 

Lj,  ;=  42. 1 29n  Distance  between  bulkheads 

R  :=  l7.3285n  Mean  radius  of  shell 
t  :=  o.263in  Thickness  of  shell 

R„  :=  R  +  -  Outer  radius  of  shell 
°  2 

Rq  =  17.46in 

Dq  :=  2  Rq  Outer  Diameter  of  shell 

Ring  Stiffeners 

:=  0.1 98n  thickness  of  web  of  ring  stiffener 
depth  :=  i.02fin  height  of  ring  stiffener 

b  :=  faying  width  of  stiffener  (from  P&S  for  I  beam  stiffener) 

bf  :=  o.763in  breadth  of  ring  stiffener 

b2 bf  - breadth  of  ring  stiffener  minus  the  web  thickness 
d  ;=  .263in  ring  Stiffener  flange  thickness 
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L|j  :=  Lj  - Unsupported  spacing  between  stiffeners 


L[,  =  2.467in 


L:=  ma>{L5,Lb) 


L  =  2.665in 


cr-=-5' 


^  j  .  2  ^  J2^ 

depth  +  b2d 
+  bjd 


first  centroidal  height  of  ring  stiffener 


C|  =  0.351in 


C2  :=  depth  -  Cj 


second  centroidal  height  of  ring  stiffener 


C2  =  0.674in 

h  :=  C]  -  d  distance  from  centroid  of  ring  stiffener  to  nearest  edge  of  flange 

Rg  :=  R  -  .5  t  -  C2  radius  to  centroid  of  ring  stiffener 
Rj  =  16.523in 

Aj  :=  depth  +  b2  d)  cross-sectional  area  of  ring  stiffener 


Aj  =  0.352in^ 

ij.  :=  Q  j(^bfc  -  bjh^  +  V^2^)  ^’^oment  of  inertia  of  ring  stiffener  about  its  centroidal  axis 
Ir=0.031in‘^ 

Rf  :=  R  -  .5  t  -  depth  Radius  to  tip  of  the  stiffener 
Rf  =  16.172in 
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Combined  Plate  and  Ring  Stiffener 


:=  mi 


Effective  length  of  shell 


Lg  =  1 .999in 


Ap:=Lt 


area  of  effective  plate 


'^sp  ■ 


area  of  plate  and  ring  stiffener 


H-  :=  depth  +  t 


height  of  combined  plate  and  ring  stiffener 


H(,  =  0.1 07ft 


B| L 


plate  length  minus  the  web  thickness 


tvv’Hc  +  Bjt  +  b2d-^2H^  — 
+  B|’t  + 


neutral  axis  of  combined  plate  and  ring  stiffener 


from  outer  fiber  of  plate  (RO) 


Cj^  =  0.4in 


:=  R  +  .5’t  -  cj^  radius  to  centroidal  axis  of  combined  ring  stiffener  and  shell 


R^  =  17.06in 


'=»  t  bf-(H,  (bf  -  ,  -  df] 


moment  of  inertia  of  combined  plate  and  shell 


Ig  =  0.187in'’ 


General  Equations 


M  =  1.248 
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1 


e  =  1.605 


Q- 


2 


Q  =  0.802 

cosh(2*Q)  “  cos(2-Q) 
sinh(2Q)  +  sin(2Q) 

N  =  0.774 


^  (sinh(Q)‘COs(Q)  +  cosh (Q)-sin(Q)) 
sinh(2'Q)  +  sin(2Q) 

G=  0.934 

sinh(2Q)  -  sin(2  Q) 
sinh(2'Q)  +  sin(2Q) 


H  =  0.41 

Inter-Stiffener  Strength  (6.19.1) 

1)  Inter-stiffener  strength  equations 
This  equates  to  axisymmetric  buckling 


A:=A. 


R  ^ 


Effective  area  of  plate  and  stiffener 


A  =  0.369in 


2 


•G 


A  +  Vt  + 


2N4L 

0 


F  =  0.267 
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yield  pressure  at  midbay  and  midplane  of  cylinder 


y  l-F 


Py  =  2039. 3psi 


von  Mises  buckling  pressure  for  a  cylinder 

This  corresponds  to  Lobar  buckling 


2.42-E- 


(,-vf. 


4  ^  X  2 

- - 0.45f— 1 

2R  I  2R } 


P^=  10369.9psi 


maximum  allowable  working  pressure  for  inter-stiffener  strength 


P  P 

m 

-  if  - <1 

2  P 

y 

f  ?  \  p 

y  m 

P^-  1 - —  if  1  <  —  <3 

^  2P  P 

V  ni; 


5  _ 

— p  if  —  >  3 

6  y  Py 


=  5.085 


P„  =  1699.4psi 


^aits  ■ 


Paits  =  1359.5psi 
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2)  longitudinal  stress 

limiting  pressure  corresponding  to  the  longitudinal  stress  at  stiffeners  reaching  yield 
No  direct  correlation  to  major  failure  modes 


A  +  + 


2-N-t-L 


2-cj  t 

?,:= - 1  + 

'  R 


Pi  =  2098psi 


P  i  :=  Pi-. 67  maximum  allowable  working  pressure  for  longitudinal  stress 


Pgl5  =  1405.69psi 


Overall  Buckling  Strength  (6.19.5) 


A2(n)  :=  n  -  1 


A,(n)  := 


A2(n)  +  ~  j- 


1  (  ^ 

'\n  “I"  A,  / 

/ 


Et  EIe'^2(") 

Pn](")  :=  ^ ; - 

^  R^L 


105 


n  >  2 

02  :=  Minimiz^Py^j  ,n) 

n2  ^  2.974  n  must  be  an  integer 

"2int'=''Ound(n2,0) 

"2int  =  3 


^2firsf  ^2int * 

^2int  “  ^  otherwise 


Pn  :=  min(pn2) 


=  4496.1psi 
^aol Pn'*^ 

P^oi  =  2.248  X  lo^psi 
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Summary 

Axisymmetric  Buckling 
Py  =  2039. ^si 
Lobar  Buckling 
10369.S|3si 

General  Instability 
Pj^  =  4496.  Ipsi 

"2int  =  ^ 
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Germanischer  Lloyd 


MPa  :=  lO^Pa 
ksi  :=  1  O^psi 

General  Definitions 


E:=3000(ksi 

modulus  of  elasticity 

Lj  :=  2.665n 

stiffener  spacing 

13:=  42.l2SIn 

length  of  cylinder  between  bulkheads  or  lines  of  support 

s  ;=  0.263in 

thickness  of  shell 

R  :=  17.3285n 

radius  to  centerline  of  shell 

V  :=  0.3 

Poisson's  Ratio 

k  :=  98.5ksi 

minimum  yield  stress  of  material 

Stiffener  Dimensions 

ef:=0,263in 

flange  thickness 

df:=  0.2825n 

width  of  flange  from  web  to  edge  of  flange 

b2  2'dj* 

d^  ;=  0.762n 

Height  of  web 

:=  O.lOSn 

web  thickness 

b:=e^ 

width  of  stiffener  ring  in  contact  with  shell 

A]  :=  ef{e^  +  2df)  +  e^-d^  cross-sectional  area  of  stiffener  ring 
A,  =  0.352in^ 
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Hstiff  =  1025in 


2  2 
1  V^stiff  +b2ef 

*  2  bjCf 

C|  =  0.351in 

•=  ^stiff  ■  ‘^1 

C2  =  0.674in 


h2:=ci  -Cf 


h2  =  0.088in 


I 


(®w  +  2  -  (2-df)  h2^  +  ^^-^2  _ 


Centroidal  Moment  of  inertia  of  ring  stiffener 


Ij  =  0.031  in'* 

e  :=  C2  +  .5  s  distance  from  stiffener  centroid  to  center  of  shell 


e  =  0.805in 

Effective  Stiffener  and  Shell 


^test  ^  +  ^J2’K^s 


^test  ^test  - 
L|  otherwise 


Effective  length  of  shell  (eqn  45a,b) 


Lg  =  2.665in 


A|e 


1  + 


12 


Lg-s 


Moment  of  Inertia  of  combined  plate  and  shell 


Ig  =  0.187in* 
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^stiffc  ^  Total  height  of  stiffener  and  plate 

^stiffc  “  l-288in 

d,:=s 

bs,:=L^-e,v 

2  2 

1  ^vv'^stiffc  ^2‘^f  r^  (^’^*stiffc  ~ 

^Ic  - - 

^  ^w'^stiffc  ^  ^2‘^f  ]  '^ 

c  I  ^  =  0.888in 

^2c  -^Hstiffc-^lc 
C2c  =  0-4  in 

Rq  :=  R  +  .5  s  -  C2c  radius  of  Stiffener  ring  centroid  including  effect  of  Lg 

Rp  =  17.06in 

r2 

A  :=  A  ] — -  Modified  area  of  stiffener  ring 

Rq 

A  =  0.363in^ 

®1  •=‘=2c  + 

2.4  Asymmetric  Buckling  (Lobar  but  not  named  this) 

7I'R 

''•1  •=“; —  eqn  19 

n  :=  2  initial  guess 


110 


Pnl(n) 


(  2  \ 
n 

- +  1 

y 


-  1  . 

12-R^-(l  -  v^) 


2 

n 


-  1  +  .5?i| 


2 


PnlW  •= 


EsPj^](n) 


R 


Given 

n  >  2 

Pg  :=  Minimiz^Pj^  j  ,n) 


n  =  10.936  n  must  be  an  integer 

o 


ngint  :=  round (pg.o) 


^gint  ”  ^  ^ 


^gfirst ' 


^gint‘^^^  *^gint  ^ 
Hgint  -  1  otherwise 


^gfirst 


=  10 


gpnme ' 


'  "gfirst  ' 

rio^ 

"gint 

11 

"gint  +  1 

^gprime 

12 

n  *  +  "i"  2 

vl3y 

V  gmt  ^ 

Pn2  ■= 


'^Pnl(^"gprimeQ^^ 

Pnl^"gprimej^ 

Pnl^"gprim^^ 


^Pnl^^^gprimeg 


Per  :=  "’'KPn2) 
p^j.  =  9269.  Ipsi 


eqn  18 


eqn  17 
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3.3.3  General  Instability 

L  :=  L3  Spacing  between  effective  stiffeners 

7t  R 

^2=  1-292 

(3„2(n)  := - - 

I  -f  0.5-X2^)-(n^-f 

E-sPn2(") 

PoW-— 5— 

(n^-  |)  e1 

P|(n) - - 

Pnla(")  •=  Po(")  +  Pl(") 

n  :=2 

Given 


n  >  2 

n2:=Minimiz^p„]3,n) 

n2  =  2.955  n  must  be  an  integer 


"2int  :=  round (n2,0) 

"2int  =  3 
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psi 


4651 

6612 

vl0240y 


2.6  Symmetric  Buckling  (axisymmetric) 

More  Definitions 

For  Elastic-Plastic  Region  (Eqns  are  in  the  program) 

Eg  ;=  E  Secant  Modulus  =  Young's  Modulus  for  elastic  region 

E(  :=  E  Tangent  Modulus  =  Young’s  Modulus  for  elastic  region 

Vp  :=  V  Poisson’s  Ratio,  elastic-plastic;  =  Poisson's  ratio  in  elastic 

shape  factor 

calculation  factor  for  symmetric  buckling 

critical  pressure,  elastic,  calculation  factor 

Past  “  8364.876psi 


jAi  -  V  J 


a  =  0.024- 


mm 


:=  a-Lj 


C^=  1.605 


Past  • 


2s^E 
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G 


P 

Past 


Ct.=  j4TTg 
Cg  :=  C5C^ 

C9  :=C5C^ 


cos^C^^sin^C^)  cosh(Cg)sinh(Cg) 


^1 

C6 

cosh(Cg^sinh(Cg) 

cos(C9)sin(C9) 

^6 

"  C, 

cosh  (Cg^sin^C^) 

sinh^Cg^cos^C^^ 

^7 

C6 

cosh  (Cg)sinh(Cg) 

cos{c^s\n^C^ 

<=6 

*  C, 

~  ^10'p2+  ''■ClO'CirF^ 

^o’(^2  ~ 

a^=  -8.275  X  10^  Pa 

*^(1)1  •=  “  ^10  *"2  “  ''■^lO'^l  r^4) 

cr^=  -1.0I3X  10^  Pa 

if  the  calculated  stress  is  >  .8  of  the  yield  stress,  then  must  use  Et,  Es,  etc 

.8  k  =  78.8ksi 


Ko^=- 


2  2 


G:  =  135.505ksi 


a  :=  -K^  g 


X 


a  =  147ksi 


a 


8i  =  4.898x  10 


a  :=  k- 


.8  +  .2tanh 


a  =  98461.1psi 
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E,  =  59.275ksi 


H3-[(l-2Vp)-(2-Vp).K„] 

H2:=[(2-Vp)-(l-2Vp).Kj 
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H2H3H4 


Pm-2931.4psi 


p  =  293  Ipsi 

iterate  on  p  until  pm  is  equal  to  p 
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Summary' 

Axisymmetric  Buckling 

Pm  =  2931.4psi 

Lobar  Buckling 
Pj.,.  =  9269.  Ipsi 
General  Instability 
pg  =  4651.4psi 
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Society  of  Naval  Architects  and  Marine  Engineers 

Submersible  Vehicle  Systems  Design 

ksi  :=  lOOC^si 

Input  Section 


Lf:=2.665n 

Length  between  frames 

Lj,  =42.125111 

Length  between  bulkheads 

R:=  17.3285n 

Radius  of  cylinder  to  centerline  of  shell 

t  :=  0.263n 

shell  thickness 

Rq  R  +  .5-t 

radius  to  outside  of  shell 

Rq=  17.46in 

E:=  30000ksi 

modulus  of  elasticity 

Gy  :=  98.Sksi 

minimum  yield  stress  of  material 

~  0.3 

Poisson's  ratio 

Ring  Stiffeners 

:=  0.1 9Sn 

thickness  of  web  of  ring  stiffener 

depth  :=  1.025n 

total  height  of  ring  stiffener 

bf  :=0.763n 

breadth  of  ring  stiffener 

d  :=  0.263n 

ring  stiffener  flange  thickness 

b2:=bf-t^ 

breadth  of  ring  stiffener  minus  the  web  thickness 
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first  centroidal  height  of  ring  stiffener 


f  2  2 

depth  +  b2d 

C  ;=.5.  - 

depth  +  b2  d 
C|  =  0.351in 

C2  :=  depth  -  c  j  Second  centroidal  height  of  ring  stiffener 

C2  =  0.674in 

h  ;=  C]  -  d  distance  from  centroid  of  ring  stiffener  to  nearest  edge  of  flange 

h  =  O.OSSin 

:=  R  -  .5 1  -  C2  radius  to  centroid  of  ring  stiffener 
R^=  I6.523in 

Aj.  ;=  (t^  depth  +  b2  d)  cross-sectional  area  of  ring  stiffener 
=  0.352in^ 

Ir-=  "  ’’2'^^  ‘w  ‘^2^)  moment  of  inertia  of  ring  stiffener  about  its  centroidal  axis 

1^=0.031  in'* 

Combined  Plate  and  Ring  Stiffener 

Ap  :=  Lf-t  area  of  effective  plate 

Asp  ■=  ^p  +  plate  and  ring  stiffener 

:=  depth  +  t  height  of  combined  plate  and  ring  stiffener 

Bj  :=  Lf  -  plate  length  minus  the  web  thickness 
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twV+Bi-t^  +  b2d(2-H,-d) 

V^c  ■'■  ■*■  ‘*■’^2 


neutral  axis  of  combined  plate  and  ring  stiffener 


from  outer  fiber  of  plate  (Ro) 


Clc  =  0-4in 


:=  R  +  .5  h  -  Cl 


radius  to  centroidal  axis  of  combined  ring  stiffener  and  shell 


=  16.972in 


-Bl  hc-‘)  +  bf-K  -  cic)  -(bf-t^)  (H^-C]^-d)- 


moment  of  inertia  of  combined  plate  and  shell 


Ig  =  0.187in 


a.  Axisymmetric  yielding 


Af  + 


B  =  0.129 


2-R  A  R 


0  =  1.605 


cosh(0)  -  cos(0) 
sinh(0)  +  sin(0) 


N  =  0.774 
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H  =  -0.833 


b.  Lobar  Buckling 


Pb  =  10369.942jDsi 


c.  General  Instability 


m:= 


71  R 


n  ;=  1 


Given 

n  >  1 

n2  :=  ,n) 

1)2  =  2.974  n  must  be  an  integer 

n2int:=round(n2,0) 

"2int  =  3 

"2first'=  "2int  ’00  if  n2int=  1 
”2int  “  ^  otherwise 

”2first  “  ^ 


11510.622^ 

4496.065 

psi 

6320.305 
9773.51  ) 


Per  :=™"(Pcr2) 
Pj,j.  =  4496.  Ipsi 
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Summary 

Axisymmetric  Buckling 
Py  =  1927.7psi 
Lobar  Buckling 
Pjj  =  10369.9psi 
General  Instability 
P(,^  =  4496.  Ipsi 
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MIT  13A  Professional  Summer  Submarine  Design  Trends 

PROGRAM  TO  COMPUTE  SUITABILITY  OF  SUBMARINE  DESIGN  PARAMETERS. 

Ref:  "Hull  Material  Trade  Off  Study",  D  Fox,  Jan  94 
Define  input  parameters: 
ksi  =  lOOQDsi 

Global  Variable  Inputs: 


e  =  O.Oin 

Eccentricity 

Material: 

ay  =  9850Q)si 

Yield  Strength 

E  =  3000Q<si 

Modulus  of  Elasticity 

V  =0.3 

Poisson’s  Ratio 

Geometry; 

34.657in 

shell  diameter 

D 

2 

shell  radius 

tf  =  .263in 

flange  tickness 

Lf  =  2.665n 

frame  spacing 

wf  =  .763in 

flange  width 

4  =  42.129in 

bulkhead  spacing 

4  =  .19Sn 

web  thickness 

tp  =  0.263n 

shell  thickness 

=  .762in 

web  height 
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Ip 

R|  :=  R  -  —  Frame  Radius 
2 

Rl  =  17.197in 

Af  :=  tf  Wf  Frame  flange,  web  area: 
bv’h\v 

A  :=  Af  +  Aw  Frame  Area  =  Flange  +  Web 

A  =  0.352in^ 

PARTI  SHELL  YIELDING 

Von  Sanden  and  Gunther  (1952) 

PNA  Section  8.4 

bv ’^p 

B  := -  Area  ratio 

A  +  tw-tp 

Slenderness  parameter 

Deflection  coefficient: 

Frame  flexability  parameter: 


P  =  1.676 


l  +  p 


Frame  deflection  parameter: 


sinh 


Bending  effect  (mem): 


V 


9 

•sinl  — 


sinh(0)  +  sin(0) 


Hm  =  -0.934 


0.5  sinh 


He:=-2- 


2 

1  -  V  y 


—  -COS  — 


cosh 


f- 


0 

sinl 


sinh(0)  +  sin(0) 


He  =  0.368 


^  sinh(0)  -  sin(0) 
sinh(0)  +  sin(0) 


Bending  efffect  near  frame: 


P  crl  •= 


fy  ^ —  Of  y 

Pq  ^  ^ 

Pj  ^  1000  psi 
5000  psi 


delta  <—  5  psi 
limit  1  psi 
conv  <—  1  psi 
j^O 

while  j  <  200 
I  for  i  e  0.,  2 


^  <W60  ^ 

^  <M)si  ^ 

^  xxso  ^ 

xxsi  ^ 

4)4)fo  ^ 


-P.R 

1 

tp 

-P-R 

1 

-P.  R 
1 

-PR 

1 

-P.R 

1 

tp 

-P.R 


[l  +  r  (H  JVI  +  v  H  eQ 
[l  +  r-i'H  M  -  v  H  Ei] 


(0.5  +  T-He' 


(0.5  -  F  He' 


_ 

r  3  11 

1  -  r- 

1  +  V- 

I  I  K 

.. 

_ 

U  -  V  ;  JJ 

■ 

^  3  f 11 

1  -  r 

1  -  V 

K 

.1-v ;  JJ 

Bending  effect  (bend) 
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r  r  /  3  1 

L  L  li-v^J  i 


-r.  .  .u.:> 

-  0  5  -  r  — = K 

L  li-v^J 

r  (  3 

- 05+ r  — - —  K 

L  li-v^J 


a, 


^  3^  ^  sy^ 

O  4  C? 


(2  2 

^SYM^\C!\  -aia2  +  a2 


f  2  2> 

^SY¥^{^3  -0304  +  04  , 


Q  2Y  +-  max 

[\^SY 

stiess^  ^  ^  SY 

test  <—  'fy  -  o  SY' 


break  if  jstressj  -  f y|  ^  limit 
Pg  <—  Pj  if  shel  dav  >  0 
P^  <-*  P^  if  shel  dav  <  0 


Pj  +-  Pg  + 

j  +"  j  +  1 


'  ^2  ^0  * 


I  out.  <-  j  ~  1 
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1994.11^ 
16  ; 


Peril  -  Pcrlp  Ipsi 


Peril  =  1994.1psi 


PART  2  LOBAR  BUCKLING 

Windenberg  Approx  of  Von  Mises  (1933) 
Assumes  n  lobes  =  Pi*D/L 
Collapse  pressure: 


PeLB  :=  “ — 

f  4 


Vd 


2.42E- 


.2,5 


PeLB=  10369.%isi 


PART  3  GENERAL  INSTABILITY 

Corrected  Bryant  Formula  (1954)  for  better  model  test  correlation 
Pressure  loading  is: 

P  :=  600Q3si 

Compute  effective  frame  spacing: 


y  =  0.717 

Compute  clear  length: 
L  :=  Lf  -  t\v 
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ti]  :=:0.5  >/T^ 
n|  =  0.266 

Web  thickness: 

n2  :=  0.5-^j  1  +  y 


t\2  =  0.655 
t^v  =  0.198in 


Fl 


0 


2  2 
cosh^nj-O)  -cos  (n2e) 

cosh  (ni'O^sin  h(ni-9)  COS  (n2-0)  ■sin(n2-0) 

n|  n2 


must  be  less  than  1 .00 


Fl  =  0.958 


F 1  is  almost  a  linear  decreasing  function  for  pressures  from  1  to  10000  psi  with  an  average  value 
of  0.91 .  This  will  be  used  in  the  following  analysis  as  the  pressure  is  the  unknown  and  therefore 
the  above  equations  cannot  be  directly  used. 

Fl  ;=  0.96 


Lgff  :=  Lc-Fi  +  t„,  Effective  shell  plate  length: 
Lgff  =  2.566in 

Theoretical  critical  lobe  number  values  are: 


i:=0..2 

f2^ 


n  := 


3 


Circumferential  Lobes 


Aeff  :=  Leff  tp  Effective  plate  area: 


Agff  =  0.675in 


2 
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m  :=  7r —  Longitudinal  Lobes: 

Ls 


m=  1.292 


I  [  I 

- ; -  -Af-  - -  -Aeff 

^  J  \  ^  J 

Aeff  +  Aw  +  Af 


Frame-plate  neutral  axis  (ref  web  centre+  toward  flange): 


yna=  -0-237in 


Uses  Parallel  Axis  Theorm:  Icor  =  I  +  Ad2 


Moments  of  inertia  for  plate,flange,web: 


^pcor  •”  ^p  ^eff 


tp  + 


2 

^wcor  Kv  +  ^w‘(yna) 


Mf  +  hvv 

Ifcor  •“  ^  “  Xna 

V  2 


Total: 


^eff  •“  ^pcor  ^wcor  ^fcor 


Moment  of  Inertia  corrected  for  neutral  axis. 


leff  =  0.185in^ 


The  critical  Elastic  General  Instability  pressure  is: 
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PcGl  :=  m'KPcGl) 


PcGl  =  4460. 2psi 

ngi:=3 

Summary 

Axisymmetric  Buckling 

Peril  =  1994.1psi 
Lobar  Buckling 
PcLB  =  10369.^si 
General  Instability 
PcGI  =  4460. 2psi 
ngi  =  3 
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Analytic  Solution 


Definitions 

ksi  :=  6.8947571  oSa 
E  :=  3000(ksi 

R:=  ]7.3285n 

D:=  2R 
Gy  98.5ksi 

L:=  2.665n 

:=  42.12gn 

II  :=  3 
t  :=  0.263in 

Ring  Stiffeners 

t^:-0.19Sn 
H:=  1.025n 

b:= 

bf:=  0.763in 

‘’2^=bf-tw 
d  :=  0.263in 
^)nash  :=  L  -  bf 


Modulus  of  Elasticity 

Radius  of  cylinder  to  centerline  of  shell 

Diameter  of  cylinder 

Yield  strength 

Length  between  stiffener  centers 
Distance  between  bulkheads 

Poison's  ratio  for  Fe/Steel 
Shell  thickness 

thickness  of  web  of  ring  stiffener 
height  of  ring  stiffener 

faying  width  of  stiffener  (from  P&S  for  I  beam  stiffener) 
breadth  of  ring  stiffener 

breadth  of  ring  stiffener  minus  the  web  thickness 

ring  stiffener  flange  thickness 
distance  from  flange  edges 
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c,:=.5 


VH^  +  b2d^ 
t^,H  +  b2d 


first  centroidal  height  of  ring  stiffener 


Cl  =  0.351in 


C2  :=  H  -  cj  second  centroidal  height  of  ring  stiffener 


C2  =  0.674in 


h:=c,  -d 


distance  from  centroid  of  ring  stiffener  to  nearest  edge  of  flange 


h  =  0.088in 


:=  R  -  .5  t  -  C2  radius  to  centroid  of  ring  stiffener 


R^=  16.523in 


:=  +  b2  d)  cross-sectional  area  of  ring  stiffener 


=  0.352in^ 


,  .  (  1  1  f.  .  3  .  .  3  3 

I^:=  -  -b2h  +VC2 

V  y 


moment  of  inertia  of  ring  stiffener  about  its  centroidal  axis 


I^=  0.03  lin^ 


effective  area  of  stiffener  eqn  [24a]  from  P&S 


A^ff  =  0.369in 


ratio  of  effective  frame  area  to  shell  area  eqn  [62]  P&S 


a  =  0.526 


P  ;=  -  ratio  of  faying  width  to  frame  spacing  eqn  [62]  P&S 


P  =  0.074 
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Combined  Plate  and  Ring  Stiffener 

Ap  :=  L  t  area  of  effective  plate 

Ajp  :=  Ap  +  A^  area  of  plate  and  ring  stiffener 

:=  H  +  t  height  of  combined  plate  and  ring  stiffener 

Hj,  =  0.107ft 

Bj  :=  L  -  plate  length  minus  the  web  thickness 

VH^^  +  B,t^+  b2-d-(2-H^-d) 

c,^  :=  .5 -  neutral  axis  of  combined  plate  and  ring  stiffener 

1  VHc  +  Brt  +  db2 

from  outer  fiber  of  plate  (RO) 

cic  =  0-4in 

:=  R  -  .5 1  -  cjg  radius  to  centroidal  axis  of  combined  ring  stiffener  and  shell 
Rj,  =  16.797in 

le  f  -  Br(^lc  -  +  bf-(H,  -  c,  -  (bf  -  -  c,,  -  d)'] 

moment  of  inertia  of  combined  plate  and  shell 
Ig  =  0.187in'* 

a.  Buckling  of  unreinforced  shells 

von  Mises  bucking  pressure: 

n  :=  1  guess  at  number  of  waves  around  the  circumference 
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Given 


n  >  1 

n  j  :=  Minimiz^Pyj.j  ,n) 

nj  =  10.906  n  must  be  an  integer 

n,in,:=  round  (npO) 

"lint= 


Pyp  =  9272. 3psi 
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c.  Asymmetric  Collapse  (Lobar  buckling) 

Windenberg  Approx  of  Von  Mises  (1933) 
Assumes  n  lobes  =  Pi*D/L 


Pci^jg  “  10369.9psi 


d.  General  Instability  of  shells  and  rings 


n  :=  1 

Given 

n  >  1 

n2  :=  Minimiz^p^Q|,n) 

n2  =  2.974  n  must  be  an  integer 
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n2int:=  round  (n2,0) 


^2int  ^ 


^2prime’ 


"2int 
"2int+  1 

"2int+  2 


Pcgi2 


PcGI^"2prim^^ 
PcGl^'’2primej^ 
PcG](^"2prim^^ 
PcGl(^"2primej^ 

Pcgibryant  •“  '^'KPcgi2) 


Pcgibryant  =4496.1psi 


Summary 

Lobar  Buckling 
PcLB“  10369.^si 
Axi  symmetric  Buckling 

Paxiy=2140.7psi 

General  Instability 
Pcgibryant  =  4496.  Ipsi 
"2int  =  ^ 
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Appendix  C:  Codes  for  Test  Cylinder  2.a 
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AMERICAN  BUREAU  OF  SHIPPING 

Rules  for  Building  and  Classing  Underwater  Vehicles,  Systems,  and  Hyperbaric  Facilities 
Definitions 


ksi  :=  lOOQjsi 


E  :=  3000Q(si 
V  :=  0.3 
Oy  :=  6550(psi 

Shell  Parameters 


Modulus  of  Elasticity 
Poisson's  Ratio 
Yield  strength 


Lj  :=  1 .366n 
:=  8.636n 

R:=  8.4179n 


t  :=  O.OSSan 


Distance  between  stiffeners 

Distance  between  bulkheads 

Mean  radius  of  shell 
Thickness  of  shell 

Outer  radius  of  shell 


Rq  =  8.461  in 

Ring  Stiffeners 

:=  0.044n 
depth  :=  0.53  ISn 

b:=V 

bf  :=  0.39Sln 
b2:=bf-t^ 


Outer  Diameter  of  shell 

thickness  of  web  of  ring  stiffener 
height  of  ring  stiffener 

faying  width  of  stiffener  (from  P&S  for  I  beam  stiffener) 
breadth  of  ring  stiffener 

breadth  of  ring  stiffener  minus  the  web  thickness 
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d  :=  0.07an 


ring  stiffener  flange  thickness 


Lj,  :=  Lg  - 


Unsupported  spacing  between  stiffeners 


Lj,  =  I.322in 


L;=  m^Lg.L^) 


L  =  1 .366in 


Cl 


2  2 
depth  +  b2*d 

tw'depth  +  b^d 


first  centroidal  height  of  ring  stiffener 


Cj  =  0.143in 


C2  ‘=  depth  -  Cj 


second  centroidal  height  of  ring  stiffener 


C2  =  0.389in 


h  Cl  -  d 


distance  from  centroid  of  ring  stiffener  to  nearest  edge  of  flange 


—  R  +  .5  t  +  C2  radius  to  centroid  of  ring  stiffener 

R^  =  8.85in 

As  :=  (t^- depth  +  b2  d)  cross-sectional  area  of  ring  stiffener 


Ag  =  0.05  lin 


MU  3  3  3 

Ir:=  -  {bf-cj  -b2h  +VC2 
V  3  y  ^ 


moment  of  inertia  of  ring  stiffener  about  its  centroidal  axis 


Ij.=  1.219X  10 


R£:=  R  +  .5t  +  depth 


Radius  to  tip  of  the  stiffener 


Rf=  8.993in 


143 


Combined  Plate  and  Ring  Stiffener 


Lg  =  1.025in 

Ap  :=  L  t  area  of  effective  plate 

^sp  ring  stiffener 

Hg  :=  depth  +  t  height  of  combined  plate  and  ring  stiffener 

Hg  =  0.05  m 

B]  :=  L  - plate  length  minus  the  web  thickness 

^  VHg^  +  +  b2  d  (2  Hg  -  d) 

‘^ic  :=  -5  J  +B,  t  +  d  b2 -  neutral  axis  of  combined  plate  and  ring  stiffener 

from  outer  fiber  of  plate  (Ro) 

Cj^  =  0.174in 

Rg  :=  R  -  .5  t  +  c,g  radius  to  centroidal  axis  of  combined  ring  stiffener  and  shell 
Rg  =  8.549in 

h  ■=  -  B,.(c,g  -  tf  +  bf.(Hg  -  c,gf  -  (bf  -  g.(Hg  -  c,g  -  df 

moment  of  inertia  of  combined  plate  and  shell 

Ig  =  O.OOSin"* 

General  Equations 

M:=-^ 

M  =  1.607 
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1 


e  =  2.066 


Q=  1.033 

cosh(2-Q)  -  cos(2-Q) 
sinh(2Q)  +  sin(2*Q) 

N  =  0.942 


Q  _  ^  (sinh(Q)-cos(Q)  +  cosh(Q)  sin(Q)) 
sinh(2Q)  +  sin(2Q) 

G- 0.835 

_  sinh(2Q)  -  sin(2>Q) 
sinh(2Q)  +  sin(2Q) 


H  =  0.631 

a  Inter-Stiffener  Strength  (6.19.1) 

1)  Inter-stiffener  strength  equations 
This  equates  to  axisymmetric  buckling 


A:=A3- 


External  stiffeners  Effective  area  of  plate  and  stiffener 


A  =  0.046in 


F:=: 


All-  -  IG 
2 


A  +  t^'t  + 


2N-t-L 

0 


F  =  0.209 
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yield  pressure  at  midbay  and  midplane  of  cylinder 


P 

y  1  _F 


Py  =  844.1psi 


This  corresponds  to  Lobar  buckling 

von  Mises  buckling  pressure  for  a  cylinder 


(,-vf  -L-o4-L 

2R  I  2R 


P  =  2947.7psi 


maximum  allowable  working  pressure  for  inter-stiffener  strength 


P  :=  if  JI?<, 

2  P 

y 

P  ^  p 

p  ,._JL 

'  I  ^  '“n.J  Py 

c  P^ 

“  P  if  —  >  3 
6  ^  P 

y 


=  3.492 


=  703.4psi 


^aits  • 


Paits  =  562.7psi 
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2)  longitudinal  stress 

limiting  pressure  corresponding  to  the  longitudinal  stress  at  stiffeners  reaching  yield 
No  direct  correlation  to  major  failure  modes 


Pj  =  848,5psi 

maximum  allowable  working  pressure  for  longitudinal  stress 


Pals  =  568.47psi 


b  Overall  Buckling  Strength  (6.19.5) 
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Given 


Summary 

Axisymmetric  Buckling 
Py  =  844.1psi 
Lobar  Buckling 
=  2947.7psi 

General  Instability 
=  6498.2psi 

”2int  =  4 
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Germanischer  Lloyd 


1988  Edition 

Chapter  2  -  Submersibles 

MPa  :=  lO^Pa 
ksi  :=  1 0^ psi 

General  Deflnitions 


E;=  SOOOOcsi 
L,  :=  1 .366n 
:=  8.636n 

s  :=  O.OSSSn 

R:=8.417gn 

V  :=  0.3 

k  :=  65.Sksi 

Stiffener  Dimensions 


Modulus  of  Elasticity 
stiffener  spacing 

length  of  cylinder  between  bulkheads  or  lines  of  support 

thickness  of  shell 
radius  to  centerline  of  shell 
Poisson’s  Ratio 

minimum  yield  stress  of  material 


ej-:=  0.078n 
df;=  .I775n 
b2  ;=  2df 
:=  0.453911 


:=  0.044n 


flange  thickness 

width  of  flange  from  web  to  edge  of  flange 

Height  of  web 
web  thickness 


b  :=  e 


W 


width  of  stiffener  ring  in  contact  with  shell 
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cross-sectional  area  of  stiffener  ring 


■•=  ®f(^w  +  2'^f)  + 

2 

A  j  =  0.051  in 
Hstiff-=  ^w+ 

H3tiff=  0.532in 

*  2  e^-H5tjff+ bjef 

Cj  =  0.143in 

^2  ■='^stiff“ 

C2  =  0.389in 

h2:=c,  -Cf 

h2  =  0.065in 

ij  :=  +  2-df)cj^  -  (2df)h2^  +  e^c2^J  Centroidal  Moment  of  inertia  of  ring  stiffener 

Ij  =  1.219X  lO^^in"* 

e  ;=  C2  +  .5  s  distance  from  stiffener  centroid  to  center  of  shell 

e  =  0.432in 

Effective  Stiffener  and  Shell 

^test  •=  b  +  yll  R-s 


3 1 .645mm 


Lg  :=  if  Lg^ggj  <  Lj  Effective  length  of  shell  (eqn  45a, b) 

L|  otherwise 
Lg  =  1.246in 
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+  ll  + 


Moment  of  Inertia  of  combined  plate  and  shell 


1^  =  7.732  X  10 


^stiffc  •  ^stiff  ® 


Total  height  of  stiffener  and  plate 


d,:=s 


bs,:=Lg-e^ 


1  ^w’^stiffc  +  ^2'^f  + 


bSi  s  (2  H5tiffc  -  s) 


®w’^stiffc  ^2’®f  bsj-s 


C|^  =  0.435in 


^2c  ’  ^stiffc  ^Ic 


C2c  =  01 83in 


Rq  :=  R  -  .5  s  +  C2c  radius  of  stiffener  ring  centroid  including  effect  of  Lf 


Rp  =  8.558in 


A  -A,- 


Modified  area  of  stiffener  ring 


A  =  0.049in^ 


®1  ■-‘=2c+ 


Cj  =  0.225in 


n  >  2 

Hg  :=  Minimiz^pj,],n) 

Pg  =  13.832  n  must  be  an  integer 

"gint  :=  round (ng,o) 

”gfirst  •”  ”gint‘^®®  ”gint  “  ^ 

Hgint  -  1  otherwise 


^gfirst 


^gprime  * 


"gfirst 
”gint 
”gint  +  ^ 
"gint  +  2 
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3.3.3  General  Instability 

L  :=  L3  Spacing  between  effective  stiffeners 

,  Tt-R 

A.2  =  3.062 

Pn2(n)  := - - 

(^n^  -  I  +  0.5  >.2^]-(n^  +  A.2^) 

E-sPn2(") 

Po(n):=-^ 

(n^  -  ll-E  lg 

Pl(n);= - - - 

Ro  L] 

Pnla(")  •=  Po^")  +  Pl(") 
n:=2 


Given 


n2  ;=  Minimiz^Pj^j3,n) 

n2  =  4.01  n  must  be  an  integer 

n2int:=round(n2,0) 

"2int  =  4 

"2int  “  ^ 

"2int 
"2int  +  ' 

"2int  +  2 

Pnla(^"2priiii«^y 

Pnla^*’2primej^ 

Pnla^'’2prime^^ 

Pnla^'’2prime^^  ^ 

pg 

Pg  =  6184.9psi 

2.6  Symmetric  Buckling  (axisymmetric) 

More  Definitions 
Default  values  for  Elastic-Plastic  Region  (Eqns  are  below) 

Eg  ;=  E  Secant  Modulus  =  Young's  Modulus  for  elastic  region 

Ej  :=  E  Tangent  Modulus  =  Young's  Modulus  for  elastic  region 

V  :=  V  Poisson’s  Ratio,  elastic-plastic;  =  Poisson's  ratio  in  elastic 
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'Ml  -V 


shape  factor 


a  =  0.06- 


C5:=aL] 


calculation  factor  for  symmetric  buckling 


Cj  =  2.066 


2s^E 

-  v^l 


critical  pressure,  elastic,  calculation  factor 


Past  =  3772.56'^si 

Must  guess  pressure  and  then  change  when  compared  to  p^  at  bottom  of  calculation 


C^  :=-  yj\  -  G 


Cy 1  +  G 


Cg CjCg 


C9 


cosh( 

%)■ 

^  -  cos  ( 

c,) 

2 

cosh  I 

I'Sinh 

(Cs 

;) 

)sin| 

cosh  (CgJ-sin^C^)  sinh(Cg)cos|c^) 

^ _ ^6 

cosh  (Cg)sinh(Cg)  cosICpjsin^Cp) 
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cos^C^^-sin^Cp)  cosh^Cg^sinh^Cg^ 


^  2  cosh^Cg^sinh^Cg^  cos^Cp^  sin^Cg) 


cosh  ^Cg)  *  sin^C^^  sinh  ^Cg^  •  cos 


F4:= 


^  2  cosh  (Cg)- sinh  ^Cg)  cos^C^)  sin(Cp) 


^10‘“ 


(  A 

1  -- 

I  2j 


s-Li 


- +  —  + 

sL|  Lj 


v'-^^r 


Cj,  ;- 


0.91 


%  := 


1  -V 

pR 


'^x’-'^oi  2  ^10’^ir^4 


0j(.= -4.675  X  10  Pa 


®(t)--  '^o'(*  ~Cio  p2+  ''■^io'^irf4) 


0^  =  -5.789x  10  Pa 


^xl'-%'  2  “‘^10‘^11'^4 


o^j  = -2.292  X  10  Pa 


^(t)i  •-  ^o'(^  “  ^io'p2  ■  ''•^lO'^l  r*^4) 


cj^j  =  -5.074X  10  Pa 
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.8k  =  52.4ksi 

If  the  calculated  stress  is  greater  than  .8k  then  must  find  Et,  Es,  etc. 


8,  =  2.799X  10 


E,  =  480.6ksi 


Vp  =  0.345 
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C|  =  0.555 


Co  =  0.163 
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7t 


trial  =  0.353 

Iterate  on  m  until  trial  <=  trial2 


m:=  1 

■  I  m  , 
trial2:=  J—(m+  1) 

trial2  =  1 


Summary 

Axisymmetric  Buckling 
p^=  1029.8psi 
Lobar  Buckling 
Pj,,.  =  2814.7psi 
General  Instability 
Pg  =  6184.9psi 
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Society  of  Naval  Architects  and  Marine  Engineers 


Submersible  Vehicle 

Systems  Design 

ksi  :=  lOOC^si 

Input  Section 

Lf:=1366n 

Length  between  frames 

;=  8.636n 

Length  between  bulkheads 

R:=  8.417911 

Radius  of  cylinder  to  centerline  of  shell 

t:=  0.085an 

shell  thickness 

Rq  R  +  .5’t 

radius  to  outside  of  shell 

Rq  =  8.461in 


E:=  30000ksi 
Gy  :=  65.Sksi 
|i  03 

Ring  Stiffeners 

:=  0.044n 

depth  :=  0.53 19n 
bf  :=0395ln 

d  ;=  0.078n 
b2  :=  bf  - 


modulus  of  elasticity 
minimum  yield  stress  of  material 

Poisson's  ratio 

thickness  of  web  of  ring  stiffener 

total  height  of  ring  stiffener 

breadth  of  ring  stiffener 

ring  stiffener  flange  thickness 

breadth  of  ring  stiffener  minus  the  web  thickness 
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first  centroidal  height  of  ring  stiffener 


2  2 
depth  +  b2d 

depth  +  b2-d 
Cj  =  0.143in 

C2  :=  depth  -  C|  second  centroidal  height  of  ring  stiffener 

C2  -  0.389in 

h  :=  Cj  ~  d  distance  from  centroid  of  ring  stiffener  to  nearest  edge  of  flange 

h  =  0.065in 

Rj.  :=  R  +  ,5*t  +  C2  radius  to  centroid  of  ring  stiffener 
R^.  =  8.85in 

:=  (t^  depth  +  b2  d)  cross-sectional  area  of  ring  stiffener 
A^  =  0.051  in^ 

-  bjh^  +  t^c2^  j  moment  of  inertia  of  ring  stiffener  about  its  centroidal  axis 

\  =  1.219X  10“^  in'* 

Combined  Plate  and  Ring  Stiffener 

Ap  :=  Lf-t  area  of  effective  plate 

Asp  •=  '^p  \  of  P^3te  and  ring  stiffener 

Hj,  :=  depth  +  t  height  of  combined  plate  and  ring  stiffener 

B,  :=  Lf  - plate  length  minus  the  web  thickness 
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neutral  axis  of  combined  plate  and  ring  stiffener 
from  outer  fiber  of  plate  (Rq) 

C]^  =  0.174in 

:=  R  -  .5  t  +  c,^  radius  to  centroidal  axis  of  combined  ring  stiffener  and  shell 
=  8.549in 

U  -  -  Br(«lc  -  f  +  “(He  -  "Ic)’  -  (“-  >w)  (He  -  =le  “  jf] 

moment  of  inertia  of  combined  plate  and  shell 
Ig  =  O.OOSin'* 


Axisymmetric  yielding 


e  =  2.066 

cosh(0)  -  cos(0) 
sinh(0)  +  sin(0) 

N  =  0.942 


twV  +  B;.t^b2d(2.H,-d) 

V^c  '^■’^2 
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Py  =  814.73^si 


Lobar  Buckling 


P(j  =  2.948x  lO^psi 

General  Instability 

tiR 
m:= - 

Lb 
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E-t 


m 


f 

2^  /  \: 

2 

1  m  f  2  2) 

n  - 

I  H - Mn  +  m  / 

V 

2  ) 

Pcrl(") 


Given 

n  >  1 
n2 

n2  =  3.957  n  must  be  an  integer 

n2int;=  round  (n2,0) 


(n^  -  i)  e  I 


"2int  =  4 


^2first  • 


n2int  l00  if  n2int=  1 
^2int  “  ^  otherwise 


^2prime' 


^  '’2first  ^ 


"2int 
"2int  +  ^ 

v"2int+2^ 


Pcr2  • 


Pcrl^”2prime^^ 

Pcrl(^"2prime,^ 


Pcrl(^"2prim^^ 
Pcrl^”2prim^^ 

Pcr:=min(pcr2) 

Pgj.  =  6498.2psi 
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Summary 

Axisymmetric  Buckling 
Py  =  814.7psi 
Lobar  Buckling 


Pj,  =  2947. 7psi 


General  Instability 
Pj.p  =  6498.^si 
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MIT  13A  Professional  Summer  Submarine  Design  Trends 

PROGRAM  TO  COMPUTE  SUITABILITY  OF  SUBMARINE  DESIGN  PARAMETERS. 
Ref:  "Hull  Material  Trade  Off  Study",  D  Fox,  Jan  94 
Define  input  parameters: 
ksi  =  100Q)si 

Global  Variable  Inputs: 

e  =  0.0-in  Eccentricity 

Material: 

Oy  =  6550Q)si  Yield  strength 

E  =  3000a<si  Modulus  of  Elasticity 

=  0-3  Poisson’s  Ratio 

Geometry: 

R  =  8.41 7gn  sliell  radius 

D  =  2  R  shell  diameter 

Lf  =  i.366n  frame  spacing 

4  =  8.636in  bulkhead  spacing 

tp  =  o.ossan  shell  thickness 

tf  =  O.OVSn  flange  tickness 

wf  =  o.39gn  flange  width 

tw  =  0.044n  web  thickness 

hw  =  0.453911  web  height 
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Compute  areas; 

tp 

Rf  :=  R  +  —  Frame  Radius 

2 

Rf  =  8.461  in 

Frame  flange,  web  area: 

Af  :=  tf- Wf 

A\v  •=  t\v’h\v 

A  :=  Af  +  Aw  Frame  Area  =  Flange  +  Web 

A  =  0.05  lin^ 


PART  1  SHELL  YIELDING 


Area  ratio 


0  :=  Lf- 


0  =  2.066 


cosh(0)  -  cos 
sinh(0)  +  sin(0) 


Slenderness  parameter; 


Deflection  coefficient: 


1  13 

^  T  ~7 - ^  Frame  flexibility  parameter: 

P  =  1 .947 
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V 

1  -- 
2 


-B 


Frame  deflection  parameter: 


r  := 


1  +  p 


r  =  0.265 


sinh 


Hm:=-2- 


Hm  =  -0.835 


He  :=  -2- 


V2 


■cos 


. ,  e)  .  fe' 

+  coshl  —  -sin  — 

^  Bending  effect  (mem): 


sinh(6)  +  sin(6) 


,0.5  sinh 


1  -V 


—  I'COS 

l2 


—  -  cosh 


\2  J  2 


sinh(0)  +  sin(6) 


Bending  effect  (bend): 


He  =  0.557 

^  sinh(6)  -  sin(0) 
sinh(0)  +  sin(0) 


Bending  efffect  near  frame: 
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Perl  ■- 


f  y  ^  a  y 

Pg  <-  1  psi 
Pj  1000  psi 


jP^  <r-  2000  psi 

delta  <-  5  psi 
limit  <—  1  psi 
conv  1  psi 
j^O 

while  j  <  20 
I  for  i  €  0..  2 


-P.  R 


_  ^  1 

[l  +  r  iHm  +  V  H  gi] 

^  (J4ko  ^ 

^P 

-P.R 

[l  +riHM-  v  Hgi] 

44bi  ^  . 

^P 

-P.  R 

<=^  xxso  ^  . 

•i0.5  +  rHE> 

-PR 


|cr  xxsi  ^ 

^  (Mfo  ^ 

^  Wifi  ^ 


^  xxfo  ^ 


^  xxfl  ^ 


‘p 


tp 

-PR 

1 

-P.R 


-P.R 

1 


•'0.5 -T  He* 


r'  3  > 

0.5  T 

1  -  r 

1  +  V 

■K 

.  2 

- 

- 

U  -  V  J 

/ 

3  \* 

D.5  T 

1  -  r- 

1 

-V 

1 

K 

,  2 

- 

\ 

,1-v  ) 

_ 

f 

^  3 

.05  ■ 

-  0.5  -  r 

-  K 

L 

<i-’ 

^  )  . 

[ 

''  3 

.0,5  - 

0.5 +  r- 

-  1  K 

2 

L 

\ 

- ;  . 

^  xxso 
^  xxsi 


a  sy 


a  1  ^  a 


o  (Jxifi 

^  xxfo 
xxfl  ; 

% 


a2<-asy^ 


^  3  ^  Of  sy^ 

CT  4  a  sy^ 
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Perl  " 

Peril  • 


t^SYM' 


{  2  2) 

-a  1-02 +  (7  2  J 


(  2  2') 

aSYF^  -o3'74  +  <^4  J 

(  fc^SYIvlV 

a  SY  <-  max 

[X^SYF  J 

stress^  <r-  o^Y 

test.  ^  ify-  asyi 
teslQ 


shel  <— 


dav  <- 


conv 
testj 

conv 


break  if  jstressj  -  fyj  <  limit 

Pq  if  shel  dav>0 

P2  ^  P|  if  shel‘dav<0 


Pj  ^  Pq  + 


^^2  ^0- 


j<- j+  1 


OUtg  4— 
1 


^  1 


conv 
out*  <-  j  -  1 


out 


^819.181'' 

^  14  y 

Pcrlplpsi 

Peril  =  819.2psi 
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PART  2  LOBAR  BUCKLING 


Windenberg  Approx  of  Von  Mises  (1933) 

Assumes  n  lobes  =  Pi*D/L 

Collapse 

pressure: 


PcLB  =  2947.7psi 

PART  3  GENERAL  INSTABILITY 

Corrected  Bryant  Formula  (1954)  for  better  model  test  correlation 
Pressure  loading  is: 

P:=  p-gDfSFgj 
P=  1786.1 3  Ipsi 

Compute  effective  frame  spacing: 


y  =  0.473 

Compute  clear  length: 
Lc  :=  Lf  -  t^y 
n  I  :=  O.S  -^  1  -  y 
ni  =0.363 
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n2  :=  0.5-^  1  +  y 
1X2  -  0-607 

2  2 
cosh(nj-0)  -  cos(n2-0) 

cosh  (ni-0)sinh(n  1-0)  cos(n2-0)sin(n2-0) 

ni  02 

F]  =  0.894  must  be  less  than  1 .00 
Fi  is  almost  a  linear  decreasing  function  for  pressures  from  1  to  2000  psi  with  an  average  value 
of  0.90.  This  will  be  used  in  the  following  analysis  as  the  pressure  is  the  unknown  and  therefore 
the  above  equations  cannot  be  directly  used. 

F]  :=  0.90 

Leff  :=  Lc'Fi  +  bv  Effective  shell  plate  length: 

Lgff  =  1 .234in 

Theoretical  critical  lobe  number  values  are:  i  :=  o..  2 
Aeff  :=  Lefftp  Effective  plate  area 

2 

Aeff  =  0.106in  I 

Circumferential  Lobes 


Longitudinal  Lobes: 
R 

m:=n — 

4 

m  =  3.062 


4 
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Ieff=7.711x  10“^  in'* 

Determine  Rf  as  radius  to  centroid  of  combined  plate  and  stiffener: 
The  critical  Elastic  General  Instability  pressure  is: 
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Min  Pressure; 


PcGl 


4^ 

2.039X  10 
8.546X  10^ 
v6.379x  10^^ 


psi 


f2^ 


n  := 


3 


PcGI  :=  min(PcGl) 


PcGI  =  6378.9psi 


Summary 

Axisymmetric  Buckling 

Peril  =  81 9.2psi 

Lobar  Buckling 
PcLB  =  2947.7psi 
General  Instability 
PcGI  =  6378. 9psi 
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Analytic  Solution 


Defintions 

ksi  :=  6.8947571  oSa 
E:=  3000(ksi 
R:=8.417gn 
D  :=  2  R 
Oy  ;=  65.5ksi 
L:=  1.36Qn 
1^  :=  8.636n 
p:=  .3 
t  ;=  0.085an 
Ring  Stiffeners 
:=  0.044n 

H:=  0.53ign 

b:=V 

bf;=  0.39gn 

b2:=bf-v 

d  :=  0.07an 
4)nash  •=  L  -  bf 


Modulus  of  Elasticity 

Radius  of  cylinder  to  centerline  of  shell 

Diameter  of  cylinder 

Yield  strength 

Length  of  supported  cylinder 
Distance  between  bulkheads 

Poison's  ratio  for  Fe/Steel 
Shell  thickness 

thickness  of  web  of  ring  stiffener 
height  of  ring  stiffener 

faying  width  of  stiffener  (from  P&S  for  I  beam  stiffener) 
breadth  of  ring  stiffener 

breadth  of  ring  stiffener  minus  the  web  thickness 
ring  stiffener  flange  thickness 
distance  from  flange  edges 
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c,:=.5. 


VH+bjd 


first  centroidal  height  of  ring  stiffener 


C|  =  0.143in 


C2  :=H-Cj 


second  centroidal  height  of  ring  stiffener 


C2  =  0.389in 


h  :=  cj  -  d  distance  from  centroid  of  ring  stiffener  to  nearest  edge  of  flange 


h  =  0.065in 


Rj.  ;=  R  +  .5  t  +  C2  radius  to  centroid  of  ring  stiffener 


Rj.  =  8.85in 


:=  (t^  H  +  b2  d)  cross-sectional  area  of  ring  stiffener 


A^=  O.OSlin^ 


moment  of  inertia  of  ring  stiffener  about  its  centroidal  axis 


1^=  1.219X  10 


^r' 


effective  area  of  stiffener  eqn  [24a]  from  P&S 


Ag^  =  0.049in 


ratio  of  effective  frame  area  to  shell  area  eqn  [62]  P&S 


a  =0.415 


ratio  of  faying  width  to  frame  spacing  eqn  [62]  P&S 


p  =  0.032 
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Combined  Plate  and  Ring  Stiffener 

Ap  ;=  L  t  area  of  effective  plate 

'^sp  1‘ing  stiffener 

Hj,  ;=  H  + 1  height  of  combined  plate  and  ring  stiffener 

Hj,  =  0.051  ft 

Bj  :=  L  - plate  length  minus  the  web  thickness 
+  B,  t^  +  b2  d  (2  H^  -  d) 

‘^ic  •=  ~H  +B,  t  +  d  b2 -  neutral  axis  of  combined  plate  and  ring  stiffener 

from  outer  fiber  of  plate  (Ro) 

Cj^  =  0.174in 

Rc  ■=  R  -  t  +  cic  radius  to  centroidal  axis  of  combined  ring  stiffener  and  shell 

Rp  =  8.549in 

:=  3  {^  <=10^  -  Bi  (c,c  -  tf  +  bf  (H^  -  c,^)^  -  (bf  -  t  J.(h^  -  c,^  -  d)^ 
moment  of  inertia  of  combined  plate  and  shell 

Ig  =  0.0 1  in^ 

a.  Buckling  of  unreinforced  shells 


guess  at  number  of  waves  around  the  circumference 
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Given 
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b.  Axisymmetric  Buckling 


Pcab  ■“ 


y  <-  0 

limits-  lOpsi 
test  <-  Opsi 
conv  ♦-  Ipsi 
j4-0 

while  j  <  20 


nr 


n2^  -  v'  +  y 


n — 


2  L 


F, 


cosh^q  I'Oj 

r  co.s( 

nypf 

cosh^ri  I'Gj 

I'Sinhj 

[n|0)  cos 

(12-6] 

l-sinl 

[12O) 

cosh 


ll  12 

(n  j-ej  sin(n2  0)  sinh^n  j  ej  cos(Ti2‘0) 


^2 


^1 


cosh^T]  j  0|  sinh(n  j  cos^H 2  ®)'S‘f’(n2‘®) 

ni  n2 

cos|n2‘®)  sin(Ti2  0j  cosher)  j-oj-sinh^Ti  j-oj 
3  ^2  ni 


>1  »- 


3 


2  cosh(n  j  Oysinh^Ti  j-G)  cos(n2’®)  s>n(n2'Oj 
H)  n2 

cosh^n  j  0j-sin|ri2-Bj  sinh^H  |■6j•cos|r] 2  8) 
"12  ni 


cosh^n  j-Gj  sinh^n  j  Oj  cos  (ij-e)  ■sin(ii2  e) 

1|  12 


('-fb 

a  +  p  +  (l-p).F, 


dcnomi  4—  A 


denom2  ^ 


F2^+F2F4(|-2m) 


0.91  1 

1  r  2  ( 

A 

[j 

.-F^J 

+  F4  A 

1  -  M  +  M  /• 

0.91 


1  -  M 


Pc2^ 


■A|F2~mF4 


0.91  ^ 

>1  1 


3 

>1  ^ 


+  denoml  -  dcnoiTi2 
break  if  |p^2  ~  |  <  limit 

Pc2 


2-E 
test  <-  p^2 
j<- j+  1 

Pc2 


out„ 


conv 


out  <- j 
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d.  General  Instability  of  shells  and  rings 


Given 

n  >  1 
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02  :=  Mininiiz^Pj,Qj,n) 

r\2  =  3.957  n  must  be  an  integer 

n2int:=  round  (n2,0) 


"2int  =  4 


^2prime' 


"2int "  * 
"2int 
"2int  +  • 
"2int  +  2 


Pcgi2 


PcGl(^"2prim^^ 
PcGl^'’2primej^ 
PcG^*’2prime2^ 
PcGl^"2prim^^ 

Pcgibryant  •“  ”’'KPcgi2) 


Pcgibryant  =  ^  10  psi 


Summary 
Lobar  Buckling 
PcLB  =  2947.7psi 
Axisymmetric  Buckling 
Paxiy= 

General  Instability 

Pcgibryant  ~  6498.2psi 

"2int  =  4 
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Appendix  D:  Codes  for  Test  Cylinder  2.c 
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AMERICAN  BUREAU  OF  SHIPPING 


Rules  for  Building  and  Classing  Underwater  Vehicles,  Systems,  and  Hyperbaric  Facilities 
Definitions 


ksi  :=  100C|5si 
E:=  3000(ksi 
V  ;=  0.3 

Oy  ;=  ISTOOQjsi 
Shell  Parameters 
Lj  :=  3.256n 
Lj.:=  115.533n 


R:=  18.882n 


t:=  0.332n 


Modulus  of  Elasticity 
Poisson's  Ratio 
Yield  strength 

Distance  between  stiffeners 
Distance  between  bulkheads 

Mean  radius  of  shell 
Thickness  of  shell 


Ro:=R  + 


t 

2 


Outer  radius  of  shell 


Rp  =  19.051in 

Do-=2Ro 


Outer  Diameter  of  shell 


Ring  Stiffeners 

:=  0.122n 


depth  :=  2.3 1  Sn 


b 


t 


vv 


bf  :=  1.552Sn 


b2  :=  bf  -  t 


w 


thickness  of  web  of  ring  stiffener 
height  of  ring  stiffener 

faying  width  of  stiffener  (from  PifeS  for  1  beam  stiffener) 
breadth  of  ring  stiffener 

breadth  of  ring  stiffener  minus  the  web  thickness 
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d  :=  0.305n 


ring  stiffener  flange  thickness 


Lf,  :=  Lg  - Unsupported  spacing  between  stiffeners 
Ljj  =  3.129in 
L:=ma)(L5,L,,) 

L  =  3.256in 

f  2  2^ 

depth  +  b2d 

Cl  :=  .5-  -  first  centroidal  height  of  ring  stiffener 

V  depth  +  bjd  j 

Cj  =  0.558in 

C2  :=  depth  -  cj  second  centroidal  height  of  ring  stiffener 

C2  =  1.757in 

h  :=  cj  -  d  distance  from  centroid  of  ring  stiffener  to  nearest  edge  of  flange 

R3  :=  R  -  .5  t  -  C2  radius  to  centroid  of  ring  stiffener 
Rg  =  16.957in 

Ag  :=  (t^  depth  +  b2  d)  cross-sectional  area  of  ring  stiffener 
Aj  =  0.729in^ 

Ij.  :=  Q  -  b2  h^  +  Vc2^  j  moment  of  inertia  of  ring  stiffener  about  its  centroidal  axis 

\  =  Q3\2m 

Rf  :=  R  -  .5  t  -  depth  Radius  to  tip  of  the  stiffener 

Rf  =  16.398in 
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Combined  Plate  and  Ring  Stiffener 


effective  plate  length 

Lg  =  2.442in 

Ap  ;=  L  t  area  of  effective  plate 

Agp  :=  Ap  +  Aj  area  of  plate  and  ring  stiffener 

Hg  ;=  depth  +  t  height  of  combined  plate  and  ring  stiffener 

Hg  =  0.221  a 

B]  :=  L  - plate  length  minus  the  web  thickness 
Kv'^c  +  B|-t  +  b2-d-^2-Hg  -  d) 

cic  •=  -5 - - — - 7- -  neutral  axis  of  combined  plate  and  ring  stiffener 

t^Hg  +  B,t  +  db2 


C|g  =  0.937in 


from  outer  fiber  of  plate  (Ro) 


Rg  :=  R  +  .5 1  -  c,g  radius  to  centroidal  axis  of  combined  ring  stiffener  and  shell 


Rg  =  18.1 14in 

•e  -  Br(‘=lc  -  *)^  +  bf-lHg  -  c,g)3  -  (bf  -  g.(Hg  -  c,g  -  df] 

moment  of  inertia  of  combined  plate  and  shell 

Ig=  1.946in'* 

General  Equations 


M  =  1.291 
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0 


■M 


0  =  1.659 


Q:= 


2 


Q  =  0.83 

cosh  (2- Q)  -  cos(2-Q) 
sinh(2-Q)  +  sin(2Q) 

N  =  0.797 


Q _  ^  (sinh(Q)  cos(Q)  +  cosh(Q)  sin(Q)) 
sinh(2-Q)  +  sin(2Q) 

G  =  0.926 

^  _  sinh(2Q)  -  sin(2*Q) 
sinh(2Q)  +  sin(2Q) 


H  =  0.435 


Inter-Stiffener  Strength  (6.19.1) 

1)  Inter-stiffener  strength  equations 

yield  pressure  at  midbay  and  midplane  of  cylinder 

This  equates  to  axisymmetric  buckling 


A:=  Aj- 


Rs 


Effective  area  of  plate  and  stiffener  (Internal  stiffeners) 


A  =  0.81  lin^ 
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A- 


G 


von  Mises  buckling  pressure  for  a  cylinder 
This  corresponds  to  Lobar  buckling 


5 


P^=  13410.8psi 


=  3509.3psi 
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maximum  allowable  working  pressure  for  inter-stiffener  strength 


^aits  • 


PaiK  =  2807.4psi 


2)  longitudinal  stress 


limiting  pressure  corresponding  to  the  longitudinal  stress  at  stiffeners  reaching  yield 
No  direct  correlation  to  major  failure  modes 


A  + 


2-N-t-L 


f  12  V 

- -  1  +  -  y-H 

R  2 

V 1  -  V  ; 


P|  =  3566.1psi 


maximum  allowable  working  pressure  for  longitudinal  stress 


Pals -Pr -67 


=  2389.29psi 


b  Overall  Buckling  Strength  (6.19.5) 


A2(n)  :=  n  -  1 
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Et  EIeA2(n) 

pni(")  ; — 

^  R^-L 

Given 
n  >  2 

t\2  :=  Minimiz^Pj^l  ,n) 

n2  =  2  n  must  be  an  integer 


n2int:=  round  (n2,0) 
”2int  “  ^ 


Pn  :=HPn2) 

Pj,  =  8642.  Ipsi 
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Summary 

Axisymmetric  Buckling 
Py  =  4211.1psi 
Lobar  Buckling 
P^=  13410.^si 
General  Instability 
Pj^  =  8642.  Ipsi 

"2int  =  2 


Germanischer  Lloyd 

1988  Edition 

Chapter  2  -  Submersibles 

MPa  :=  loSa 


ksi  :=  10 


psi 


General  Deflnitions 


E:=  SOOOQcsi 
Lj  :=3.256n 


Modulus  of  Elasticity 
stiffener  spacing 


1^:=  1  15.532n 


s  :=  0.337in 


R:=  18.882n 


V  :=  0.3 
k:=  I57ksi 


length  of  cylinder  between  bulkheads  or  lines  of  support 

thickness  of  shell 
radius  to  centerline  of  shell 
Poisson's  Ratio 

minimum  yield  stress  of  material 


Stiffener  Dimensions 


e^:=  0.305n 
df:=  0.7125n 

:=  2.01in 

:=  0.121n 


flange  thickness 

width  of  flange  from  web  to  edge  of  flange 

Height  of  web 
web  thickness 

width  of  stiffener  ring  in  contact  with  shell 
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A  j  :=  Cf  (e^  +  2  df)  +  cross-sectional  area  of  stiffener  ring 

Aj=0,729in^ 

Hstiff=2-315in 


1  ^w'^stiff  +b2ef 


*  2  e^H5tjff+ bjCf 


C|  =  0.558in 


^2  •-  ^stiff  “ 


C2  =  1.757in 


hi Cj  -  Cf 


h2  =  0.253in 


h  ■■= 


+  2df)cj^  -  (2df)h2^  +  e^c2^J  Centroidal  Moment  of  inertia  of  ring  stiffener 


Ij  =0.312in 


e  :=  C2  +  .5  s  distance  from  stiffener  centroid  to  center  of  shell 


e  =  1 .925in 

Effective  Stiffener  and  Shell 

Latest  :=  b  + 


Letest  =  93.838iiim 


^test  ^test  - 


Effective  length  of  shell 


L|  otherwise 


Lg  =  3.256in 
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lg=  1.946in 


Moment  of  Inertia  of  combined  plate  and  shell 


^stiffc  •  ^stiff  ® 


Total  height  of  stiffener  and  plate 


Hstiffc  =  2.652in 


d,:=s 


bs,;=Le-e^ 


2  2 
2  +  b2'ej'  +  bsys 


C|^  =  l,715in 


^2c '  ^stiffc  ^Ic 


C2c  =  0.937in 


Rq  :=  R  +  .5*s  -  C2c  radius  of  stiffener  ring  centroid  including  effect  of 


Rq  =  18.1 14in 


A:=A|.- 


Modified  area  of  stiffener  ring 


A  =  0.792in^ 


ej  ;=  C2c  +  .5  s 


Cj  =  l.lOSin 
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2.4  Asymmetric  Buckling  (Lobar  but  not  named  this) 


n  :=  2  initial  guess 

-  1  + 

12  R^-(]  -  v^) 

■  '  2,2 
n  -  1  +  .5>.j 

EsPn](n) 

Pnl(")  := - - 

Given 


n  >  2 

Hg  :=Mininiiz^p^],n) 

Hg  =  1 0. 1 7 1  n  must  be  an  integer 

"gint  :=  round (ng,0) 


"gint 


10 


"gfirst  • 


"gint 

"gint 


■100  if  ngi„t=  1 

-  1  Otherwise 


^gprime 


^gfirst 
"gint 
”gint  ^ 
"gint  +  2 


^Pnl(^"gprimeQj 

Pnl^^gprimcj  j 

Pn2  •=  /  N 

Pnl(^"gprime2j 


Pcr:=  min(pn2) 


p„  =  12139.3psi 


3.3.3  General  Instability 


L:=L3 


Spacing  between  effective  stiffeners 


Pn2(")  := 


n  -  1  +  0.5- ^2  J'l  ^2 


2U  2  .2 


Po(n)  := 


EsPn2(") 


Pl(n)  := 


(n^  -  i)  e  I 


\  Li 


Pnla(")  •=  Po(")  +  Pl(") 


n2:=Minimiz^p^,3,n) 


n2  =  2  n  must  be  an  integer 
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1.768X  10^ 
9.702X  10^ 
2.41 8x 
4.525  X  lo'* 


psi 


2.6  Symmetric  Buckling  (axisymmetric) 

More  Definitions 

For  Elastic-Plastic  Region  (Eqns  are  in  the  program) 

Eg  :=  E  Secant  Modulus  =  Young’s  Modulus  for  elastic  region 

Ej  :=  E  Tangent  Modulus  =  Young's  Modulus  for  elastic  region 

:=  V  Poisson’s  Ratio,  elastic-plastic;  =  Poisson's  ratio  in  elastic 

shape  factor 

a  =  0.02— 
mm 
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198 


cosh  (Cg)  -  sin  sinh  ^Cg)  ■  cos 


¥4:= 


C-T 


^  2  cosh^Cg^sinh^Cgj  cos^C^ysin^C^) 


^lo'-' 


1  -  - 

2y 


sLj 


A  b  /", 

+  —  +  1  ] - 

My 


sL]  Lj 


•Fl 


Cj]  :- 


0.91 


1  -V 


-p-R 


‘^x--  %i  2  ^io'^ir^4^ 


2“Mo'^ir*'4 


c^=  -lJ68x  10  Pa 


cr^j^  =  -1.259x  10  Pa 


Ko:=- 


2  2 
+^(1) 


a;=  176.383ksi 


c7;=-Koax 


a  =  182.6ksi 
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E^  =  1.067X  lo'^^Pa 
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trial  := 


“i-Li 

7r 

trial  =  0.413 

iterate  on  m  until  trial  <=  trial2 


m:=  1 


rial2:= 


m 


trial2:=  f  — *(m+  1) 
2 


trial2=  1 


Pm*  Paa*^0* 


V  H-m  J 


1  f  7t-m  ^ 


p  =  4566.9psi 


p  =  456'^si 


iterate  on  p  until  equals  p 


Summary 

Axisymmetric  Buckling 
Pjj^  =  4566.9psi 
Lobar  Buckling 
=  I2l39.3psi 
General  Instability 
Pg  =  9702.4psi 
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Society  of  Naval  Architects  and  Marine  Engineers 

Submersible  Vehicle  Systems  Design 

ksi  :=  1  OOQ^si 

Input  Section 

Lf:=3.256n 

Length  between  frames 

Lj,:=  115.53an 

Length  between  bulkheads 

R:=  18.883n 

Radius  of  cylinder  to  centerline  of  shell 

t  :=  0.33'^n 

shell  thickness 

Rq  !=  R  +  .5*t 

radius  to  outside  of  shell 

Rq=  19.051in 

E:-  30000ksi 

modulus  of  elasticity 

(jy  :=  157  ksi 

minimum  yield  stress  of  material 

:=  0.3 

Poisson's  ratio 

Ring  Stiffeners 

:=  0.12'^n 

thickness  of  web  of  ring  stiffener 

depth  :=2.315n 

total  height  of  ring  stiffener 

bf  :=  1.552n 

breadth  of  ring  stiffener 

d  :=  0.305n 

ring  stiffener  flange  thickness 

breadth  of  ring  stiffener  minus  the  web  thickness 
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first  centroidal  height  of  ring  stiffener 


^  2  2 
depth  +  b2d 

c,:=.5-  - 

t^depth  +  bjd 

Cj  =  0.558in 

C2  :=  depth  -  c  j  second  centroidal  height  of  ring  stiffener 

C2=  1.757in 

h  :=  cj  -  d  distance  from  centroid  of  ring  stiffener  to  nearest  edge  of  flange 

h  =  0.253in 

:=  R  -  .5 1  -  C2  radius  to  centroid  of  ring  stiffener 

R^=  16.957in 

Aj.  :=  (t^y- depth  +  bjd)  cross-sectional  area  of  ring  stiffener 

=  0.729in^ 

3  3  3^ 

bf  cj  -  b2  h  +  t^y  C2  J  moment  of  inertia  of  ring  stiffener  about  its  centroidal  axis 

lr=0.312in'‘ 

Combined  Plate  and  Ring  Stiffener 

Ap  :=  Lf-t  area  of  effective  plate 

Asp  •=  Ap  +  stiffener 

Hj,  :=  depth  +  t  height  of  combined  plate  and  ring  stiffener 

B]  ;=  Lf  -  plate  length  minus  the  web  thickness 


204 


neutral  axis  of  combined  plate  and  ring  stiffener 
from  outer  fiber  of  plate  (Ro) 

Cjj,  =  0.937in 

Rp  :=  R  +  .5t  -  cjj,  radius  to  centroidal  axis  of  combined  ring  stiffener  and  shell 

Rj.  =  18.1 14in 

U  ■  'w)'^  -  «1C  -  ■if] 

moment  of  inertia  of  combined  plate  and  shell 
Ig  =  1 .946in'^ 


a.  Axisymmetric  yielding 


e  =  1.659 

cosh(9)  -  cos{6) 
sinh(0)  +  sin(0) 

N  =  0.797 


yHg  +B]-t  +  b2  d  (2-Hg  -  d) 
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Py  =  3863.5  Ipsi 


b.  lobar  buckling 


P|j  =  1.341x  lO^psi 


c.  General  Instability 

JiR 
m:= - 

Lb 

n  :=  1 


206 


207 


Summary 

Axisymmetric  Buckling 
Py  =  3863.5psi 
Lobar  Buckling 


Pj,  =  13410.^si 


General  Instability 
=  8642. Ipsi 
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MIT  13A  Professional  Summer  Submarine  Design  Trends 

PROGRAM  TO  COMPUTE  SUITABILITY  OF  SUBMARINE  DESIGN  PARAMETERS. 
Ref:  "Hull  Material  Trade  Off  Study",  D  Fox,  Jan  94 
Define  input  parameters: 

ksi  =  100Q)si 

Global  Variable  Inputs: 


e  =  0.0  in 

Eccentricity 

Material: 

Gy  =  15700Q)si 

Yield  Strength 

E  =  3000Cksi 

Modulus  of  Elasiticity 

d 

III 

Poisson’s  Ratio 

Geometry: 

R=  18.882n 

shell  radius 

D-2R 

shell  diameter 

Lf  =  3.256n 

frame  spacing 

4  =  115.533n 

bulkhead  spacing 

tp  =  0.332n 

shell  thickness 

tf  =  0.305n 

flange  tickness 

wf  =  1.552n 

flange  width 

4  =  0.1 22n 

web  thickness 

-  2.01in 

web  height 
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Compute  areas: 

tp 

Ri  :=  R  -  —  Frame  Radius 

2 

R]  =  18.713in 

Frame  flange,  web  area: 

Af  tf- Wf 

Aw  *= 

A  :=  Af  +  Aw  Frame  Area  =  Flange  +  Web 

A  =  0.729in^ 


PARTI  SHELL  YIELDING 

bv'tp 
A  +  tw'tp 

B  =  0.055 


cosh(9)  -  005(9) 
sinh(9)  +  sin(9) 


Slenderness  parameter: 


Deflection  coefficient: 


2-N 

A  +  tw'tp 


Frame  flexibility  parameter: 


p  =  1.366 
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PART  3  GENERAL  INSTABILITY 

Corrected  Bryant  Formula  (1954)  for  better  model  test  correlation 
Pressure  loading  is: 

P  :=  p-g-DfSFgi 

P=  1.674X  lO^psi 

Compute  effective  frame  spacing: 


7  =  0.145 
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Compute  clear  length: 

:=  Lf  ~  hv 
n\  :=  0.5-^  1  -  y 
n|  =0.462 
Web  thickness: 
r\2  := 

112  =  ^>-535 
t^v  =  0.127in 

2  2 
cosh^n|-0^  -  cos  ("20) 

cosh(n|-9)sinh(ni-6)  cos  ("2-0)  sin(n2-0) 

"I  "2 

Fi  =  0.959  must  be  less  than  1 .00 

FI  is  almost  a  linear  decreasing  function  for  pressures  from  1  to  10000  psi  with  an  average  value 
of  0.96.  This  will  be  used  in  the  following  analysis  as  the  pressure  is  the  unknown  and  therefore 
the  above  equations  cannot  be  directly  used. 

F]  :=0.96 

Effective  shell  plate  length: 

Leff:=  Lc-F]  +  bv 

Leff=3.131in 

Theoretical  critical  lobe  number  values  are;  i  :=  o.,  2 
Aeff  :=  Leff  ‘p  Effective  plate  area: 

Aeff=  1.055in^ 
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Circumferential  Lobes 
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Total: 


leiT  :=  Ipcor  +  Iwcor  +  Ifcor  Moment  of  Inertia  corrected  for  neutral  axis. 


leff=  1.92in'* 

The  critical  Elastic  General  Instability  pressure  is: 


Min  Pressure: 


f  3^ 

8.536X  10 

PcGI  = 

2.1 07x  lO'^ 

psi  n  := 

3 

^3.942x  lo"*; 

^4, 

PcGI  ••=  min(Pj.Gi) 

^cGI  =  8535. 8psi 

Hgi  :=  2 

Summary 

Axisymmetric  Buckling 
Peril  =  3712psi 
Lobar  Buckling 
PcLB=  I3410.^si 
General  Instability 
PcGI  =  8535. 8psi 

ngi=  2 
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Analytic  Solution 


Definitions 

ksi  ;=  6.8947571  oSa 
E  :=  3000Qcsi 

R:=  18,882n 

D  :=  2  R 
Gy  :=  157ksi 
L:=3.256n 

:=  115.532n 

^:=  3 
t:=0.332n 

Ring  Stiffeners 

:=  0.1 22n 

H:=2.315n 

b:= 

bf  :=  1.552n 
b2  :=  bf  - 
d  :=  0.305n 
Mjnash  :=  L  -  bf 


Modulus  of  Elasticity 

Radius  of  cylinder  to  centerline  of  shell 

Diameter  of  cylinder 

Yield  strength 

Length  of  supported  cylinder 
Distance  between  bulkheads 

Poison's  ratio  for  Fe/Steel 
Shell  thickness 

thickness  of  web  of  ring  stiffener 
height  of  ring  stiffener 

faying  width  of  stiffener  (from  P&S  for  I  beam  stiffener) 
breadth  of  ring  stiffener 

breadth  of  ring  stiffener  minus  the  web  thickness 

ring  stiffener  flange  thickness 
distance  from  flange  edges 
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c,:=  .5- 


VH^  + 


first  centroidal  height  of  ring  stiffener 


C|  =  0.558in 


C2  :=  H  -  C|  second  centroidal  height  of  ring  stiffener 


C2  =  1.757in 


h  :=  Cl  ~  d 


distance  from  centroid  of  ring  stiffener  to  nearest  edge  of  flange 


h  =  0.253in 


Rj.  :=  R  -  .5  t  -  C2  radius  to  centroid  of  ring  stiffener 


R^=  16.957in 


Aj.  :=  ( VH  +  b2  d)  cross-sectional  area  of  ring  stiffener 


=  0.729in^ 


^  "  f  3]  -  ^2  ^^  +  VC2^)  moment  of  inertia  of  ring  stiffener  about  its  centroidal  axis 


Ij=0.3I2in'‘ 


^eflf  • 


effective  area  of  stiffener  eqn  [24a]  from  P&S 


AgiY  =  0.81  lin 


ratio  of  effective  frame  area  to  shell  area  eqn  [62]  P&S 


a  =  0.739 


P  -  ratio  of  faying  width  to  frame  spacing  eqn  [62]  P&S 


P  =  0.039 
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Combined  Plate  and  Ring  Stiffener 


Ap  :=  L  t  area  of  effective  plate 

A  :=  A  +  A.  area  of  plate  and  ring  stiffener 

sp  p  1 

Hj,  ;=  H  +  t  height  of  combined  plate  and  ring  stiffener 

=  0.221ft 

Bj  ;=  L  - 1.^  plate  length  minus  the  web  thickness 


C]c:=  .51 


VH/  +  Blt"  +  b2-d.(2.H,-d) 


V^c  +  Bl’t  +  ‘^^2 


neutral  axis  of  combined  plate  and  ring  stiffener 
from  outer  fiber  of  plate  (Ro) 


=  0.937in 


R^:=R  +  .St-cjc 


radius  to  centroidal  axis  of  combined  ring  stiffener  and  shell 


Rj,=  18.1 14in 

le  -  (bf  -  g-(H,  -  Cj,  -  d)'] 

moment  of  inertia  of  combined  plate  and  shell 


Ig=  1.946in‘^ 


a.  Buckling  of  unreinforced  shells 

von  Mises  bucking  pressure: 

n  :=  1  guess  at  number  of  waves  around  the  circumference 
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Purl(")  := 


Given 


n  >  1 

n|  :=  Minimiz^Pyj.| 

=  10.14  n  must  be  an  integer 

n,i„t:=  round  (n,.0) 


4 

Pyj.  =  1.214x  10  psi 
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^40803^ 

Pcab  “  - 

\  2  y 

Paxiy-=  Pcabfj'P^' 

Paxiy  “  40803psi 

c.  Asymmetric  Collapse  (Lobar  buckling) 

Windenberg  Approx  of  Von  Mises  (1933) 
Assumes  n  lobes  =  Pi*D/L 
Collapse  pressure: 


Given 

n  >  1 
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n2  :=  Minimiz^p^Qj.n) 

n2  =  1.654  n  must  be  an  integer 


n2int:=round(n2,0) 


™2int  ^ 


^2prime‘ 


"2int 
"2int  +  * 
V"2int+2; 


Pc2i2  • 


PcG](^"2prim^j 

PcGl(^"2primej^ 
PcGl(^"2prim^j 
PcGl^'^2prim^^ 

Pcgibryant  •=  "’'KPcgi2) 


Pcgibryant  =  8-642x  10  psi 


Summary 

Lobar  Buckling 
PcLB  ~  13410.%)si 
Axisymmetric  Buckling 

Paxiy=4080.3psi 

General  Instability 

Pcgibryant  =  8642.1psi 

"2int  =  2 
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